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1 Introduction 

A major unsolved problem in theoretical physics is to reconcile the classi- 
cal theory of general relativity with quantum mechanics. These lectures will 
deal with an attempt to describe quantum gravity as a path integral over 
geometries. Such an approach has to be non-perturbative since gravity is a 
non-renormalizable quantum field theory when the dimension of spacetime is 
four. In that case the dimension of the gravitational coupling constant G is 
-2 in units where H = 1 and c = 1 and the dimension of mass is one. Thus 
conventional, perturbative quantum field theory is only expected to be good 
for energies 

< 1/G. (1) 

That is still perfectly good in all experimental situations we can imagine in 
the laboratory, but an indication that something "new" has to happen at 
sufficiently large energy, or equivalently, at sufficiently short distances. It is 
possible, or maybe even likely, that a breakdown of perturbation theory when 
(1) is not satisfied indicates that new degrees of freedom should be present in 
a theory valid at higher energies. Indeed, we have a well known example in 
the electroweak theory. Originally the electroweak theory was described by a 
four-fermion interaction. Such a theory is not renormalizable and perturbation 
theory breaks down at sufficiently high energy. In fact it breaks down unless 
the energy satisfies (1) with the gravitational coupling constant G replaced 
the coupling constant Gf in front of the four- Fermi interaction (since Gp 
also has mass dimension -2). The breakdown reflects the appearance of new 
degrees of freedom, the W and the Z particles, and the four-Fermi interaction 
is now just an approximation to the process where a fermion interacts via W 
and Z particles with another fermion. The corresponding electroweak theory 
is renormalizable. 

When it comes to gravity there seems to be no "simple" fix like the one 
just described. However, string theory is an example of a theory which tries 
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to solve the problem by adding (infinitely many) new degrees of freedom. 
Loop quantum gravity is another approach to quantum gravity which tries 
to circumvent the problem of non-renormalizability by introducing rules of 
quantization which are unconventional from a perturbative point of view. The 
point of view taken here in these lectures is much more mundane. In a sum- 
over-histories approach we will attempt to define a non-perturbative quantum 
field theory which has as its infrared limit ordinary classical general relativity 
and at the same time has a nontrivial ultraviolet limit. From this point of 
view it is in the spirit of the renormalization group approach, first advocated 
long ago by Weinberg [1], and more recently substantiated by several groups 
of researchers [2] . 

To understand the possibility of a nontrivial ultraviolet fixed point let us 
first apply ordinary perturbation theory to quantum gravity in a regime where 

(1) is satisfied. One can in a reliable way calculate the lowest-order quantum 
correction to the gravitational potential of a point particle: 

G G{r) ^, . G \ 167 

G{r) = G\l-uj— + --- ] , u; = —. 2 
r r V ^ / 307r 

Thus the gravitational coupling constant becomes scale-dependent and trans- 
ferring from distance to energy we have 

GiE)=Gil-.GE^ + :.)^^^. (3) 

It should be stressed that the scenario described in (2) and (3) is completely 
standard in quantum field theory. Let us take the simplest quantum field the- 
ory relevant in nature: quantum electrodynamics. The electron as we observe 
it in low-energy scattering experiments is screened by vacuum polarization: 
virtual electron-positron pairs, created out of the vacuum and annihilated 
again so fast that one has consistency with the energy-time uncertainty rela- 
tions, act like dipolcs and the observed charge becomes less than the "bare" 
charge. To lowest order (one loop) (2) and (3) are replaced by 

e'^(r) ( \ 

> — e(r)=e 1-— In(mr) H ), mr«;l. (4) 

r r \ on'' / 



e\E) = 1 + — ln(£/m) +•••)« , (5) 

The last expression in eq. (5) is precisely the renormalization group-improved 
formula for the running coupling constant in QED. In fact, it is easy to calcu- 
late the so-called (one-loop) /3-function for QED from the first equation in (5) 
and use this /3- function to obtain the expression on the left-hand side of (5). 
Contrary to (3) it breaks down for sufficiently high energy: it has a so-called 
Landau pole and it reflects that we expect the interactions of QED to be in- 
finitely strong at short distances. We do not believe that such quantum field 
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theories really exist as "stand-alone" theories. They either require an explicit 
cut-off, which will then enter in the obscrvablcs at high energies, or they have 
to be embedded in a larger theory without a Landau pole. 

Assume that the last expression in (3) was exact for all E (which it is not). 
One can argue that one should use G{E) in eq. (1) rather than G, in which 
case one obtains 

Thus, assuming (3) it suddenly seems as if quantum gravity had become a 
reliable quantum theory at all energy scales, the reason being that the effective 
coupling constant G{E) becomes weaker at high energies. The behaviour (3) 
can be described in terms of a /3-function for quantum gravity. Introduce the 
dimensionless coupling constant G{E): 

G{E) = G{E)E^. (7) 

From (3) it follows that G{E) satisfies the following equation: 

E^=f3{G), 0{G) = 2G-2uG'. (8) 

The zeros of /3(G) determine the fixed points of the running coupling constant 
G{E). The zero at G = is an infrared fixed point: for E ^ the coupling 
constant G{E) and correspondingly the coupling constant G{E) G. 
The zero at G = l/w is an ultraviolet fixed point : for E —* oo the coupling 
constant G{E) 1/uj (and GiE) ^ as 1/{ujE^)). 

While there is no compelling reason to take the above arguments very 
seriously, since they are based on a one-loop calculation, the claim in [2] is 
that a more careful analysis using the so-called exact renormalization group 
equations or a systematic 2 + e expansion confirms the picture of a non-trivial 
UV fixed point. A typical result of the renormalization group calculation is 
shown in Fig. 1 in the case where the effective action has been "truncated" to 
just include the Einstein term and a cosmological term. The figure has a UV 
fixed point from where the flow to low energies starts and an infrared fixed 
point (the origin in the coupling constant coordinate system). 

Where do we meet such scenarios (IR fixed points at zero coupling and a 
nontrivial UV fixed point for a suitably defined dimensionless coupling con- 
stant)? Assume one has an asymptotically free field theory in d dimensions, 
i.e. g( = is an UV fixed point. The strong interactions in four-dimensional 
fiat spacetime, quantum chromo dynamics (QCD), are such a theory. In high- 
energy scattering experiments the effective, running coupling constant goes 
to zero. This has been beautifully verified in high-energy experiments. The 
non-linear sigma model in two-dimensional spacetime is another model. It 
plays a very important role in string theory, but even before that was ex- 
tensively studied as a toy model of QCD since it is asymptotically free (i.e. 
the running, effective coupling constant goes to zero at high energies). The 
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Fig. 1. The flow of dimensionless gravitational and cosmological coupling constants 
G and A from a UV fixed point. One line flows to the infrared Gaussian flxed point 
where G and A are zero. 

asymptotically free theories have a negative /3-function. This is what makes 
the running coupling constant go to zero at high energies. Change now (arti- 
ficially) the dimension of spacctimc infinitcsimally from dtod + e. Then the 
/3-function to lowest order in e will change as follows; 

(id{9)^ (id+e^eg + (3g{g), (9) 

and the situation is as shown in Fig. 2:^ = changes from an UV fixed point to 
an infrared fixed point while the new UV fixed point will be displaced to finite 
positive value gd^) of g, a value which goes to zero when e goes to zero. In 
the case of gravity we have formally a renormalizable theory when d = 2, the 
dimension where the gravitational coupling constant G is dimensionless. One 
can show that the two-dimensional theory can be viewed as asymptotically free 
(see eq. (113) below) To apply eq. (9) to four-dimensional quantum gravity 
starting with a renormalizable theory of gravity means that e has to be two, 
which is not very small. Thus the considerations above make little sense at a 
quantitative level, but the use of the exact renormalization group indicates, 
as mentioned, that the qualitative picture is correct. 

The discussion above shows that there might be a chance that quantum 
gravity can be defined as an ordinary quantum field theory at a non-trivial 
ultraviolet fixed point. It clearly requires non-perturbative tools to address the 
question of the existence of such a fixed point and to analyze the properties of 
the field theory defined by approaching the fixed point. One way to proceed is 
by using a lattice regularization of the quantum field theory in question. The 
lattice provides an UV regularization of the quantum field theory, namely, the 
inverse lattice spacing 1/a. The task is then to define a suitable "continuum" 
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Fig. 2. The change in the beta function f5{g) in an asymptotically free theory when 
the dimension changes from the critical dimension d, where the coupling constant g 
is dimensionless, to d + e. 



limit of this lattice theory. The procedure used is typically as follows: let 
0{xn) be an observable, Xn denoting a lattice point. We write Xn = an, n 
measuring the position in integer lattice spacings. One can then obtain, either 
by computer simulations or by analytical calculations, the correlation length 
^((?o) in lattice units, from 

- log(0(.T„)0(|/„)) - |n - m|/C(5o) + o(|n - m|). (10) 

A continuum limit of the lattice theory may then exist if it is possible to fine- 
tune the bare coupling constant go of the theory to a critical value such 
that the correlation length goes to infinity, ^{go) — > oo. Knowing how C(ffo) 
diverges for go gQ determines how the lattice spacing a should be taken to 
zero as a function of the coupling constants, namely, 

^(fo) oc -. a{go) oc \go - gX- (H) 

150 — 5ol 

This particular scaling of the lattice spacing ensures that one can define a 
physical mass niph by 

mpha{go) = l/C(ffo) 

such that the correlator {0{xn)0{yrn)) falls off exponentially like e~'"p'»l^"~''™l 
for go ^ .9o when — 7y,„|, but not |n — m|, is kept fixed in the limit go g^. 
Thus we have created a picture where the underlying lattice spacing goes to 
zero while the physical mass (or the correlation length measured in physical 
length units, not in lattice spacings) is kept fixed. This is the standard Wilso- 
nian scenario for obtaining the continuum (Euclidean) quantum field theory 
associated with the critical point of a second order phase transition (for 
second-order phase transitions there exists a correlation length which diverges, 
usually associated with the order parameter characterizing the transition). 

We would like to apply a similar approach to quantum gravity, and thus 
obtain a new way to investigate if quantum gravity can be defined non- 
perturbatively as a quantum field theory. The predictions from such a theory 
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could then be compared with the renormalization group predictions related 
to the asymptotic safety picture described above. It should be mentioned that 
the asymptotic safety picture is not the only suggestion for a continuum quan- 
tum theory of gravity using only "conventional" ideas of quantum field theory. 
Very recently two other scenarios have been suggested. One is called Lifshitz 
gravity [3] and is a theory where the non-renormalizability of the Einstein- 
Hilbert theory is cured by adding higher-order spatial derivatives in a way 
somewhat similar to what Lifshitz did many years ago in statistical models. 
In fact, the setup of the theory has some resemblance with the lattice-theory 
setup of "Causal Dynamical Triangulations (CDT)" , to be described below, 
since a time foliation is assumed and the infrared limit is that of GR. How- 
ever, contrary to Lifshitz gravity, we do not attempt to put in higher spatial 
derivatives in the lattice theory. However, when a continuum limit in the 
lattice theory is taken in a specific way which is not entirely symmetric in 
space and time one cannot rule out that higher spatial derivatives can play 
a role. The other model goes by the name of "scale-invariant gravity" [5, 4]. 
It modifies gravity into a renormalizable theory by introducing a scalar de- 
gree of freedom in addition to the transverse gravitational degrees of freedom. 
Also this model has interesting features not incompatible with the results of 
computer simulations using the CDT lattice model. 

As already mentioned, we will use a lattice approach known as causal 
dynamical triangulations (CDT) as a regularization. In Sec. 2 we give a short 
description of the formalism, providing the definitions which arc needed later 
to describe the measurements. CDT establishes a non-perturbative way of 
performing the sum over four-geometries (for more extensive definitions, see 
[6, 8]). It sums over the class of pieccwise linear four-gcomctries which can 
be assembled from four-dimensional simplicial building blocks of link length 
a, such that only causal spacetime histories are included. The challenge when 
searching for a fi,eld theory of quantum gravity is to find a theory which behaves 
as described above, i.e. as in eq. (11). The challenge is three-fold: (i) to find 
a suitable non-perturbative formulation of such a theory which satisfies a 
minimum of reasonable requirements, (ii) to find observablcs which can be 
used to test relations like (10), and (iii) to show that one can adjust the 
coupling constants of the theory such that (11) is satisfied. Although we will 
focus on (i) in what follows, let us immediately mention that (ii) is notoriously 
difficult in a theory of quantum gravity, where one is faced with a number 
of questions originating in the dynamical nature of geometry. What is the 
meaning of distance when integrating over all geometries? How do we attach 
a meaning to local spacetime points like Xn and y„? How can we define at all 
local, diffeomorphism-invariant quantities in the continuum which can then 
be translated to the regularized (lattice) theory? - What we want to point out 
here is that although (i)-(iii) are standard requirements when relating critical 
phenomena and (Euclidean) quantum field theory, gravity is special and may 
require a reformulation of (part of) the standard scenario sketched above. We 
will return to this issue later. 
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Our proposed non-perturbative formulation of four-dimensional quantum 

gravity has a number of nice properties. 

First, it sums over a class of piecewise linear geometries. The characteristic 
feature of piecewise linear geometries is that they admit a description without 
the use of coordinate systems. In this way we perform the sum over geometries 
directly, avoiding the cumbersome procedure of first introducing a coordinate 
system and then getting rid of the ensuing gauge redundancy, as one has to 
do in a continuum calculation. Our underlying assumptions are that 1) the 
class of piecewise linear geometries is in a suitable sense dense in the set of all 
geometries relevant for the path integral (probably a fairly mild assumption), 
and 2) that we arc using a correct measure on the set of geometries. This is 
a more questionable assumption since we do not even know whether such a 
measure exists. Here one has to take a pragmatic attitude in order to make 
progress. We will simply examine the outcome of our construction and try to 
judge whether it is promising. 

Secondly, our scheme is background-independent. No distinguished geom- 
etry, accompanied by quantum fluctuations, is put in by hand. If the CDT- 
regularized theory is to be taken seriously as a potential theory of quantum 
gravity, there has to be a region in the space spanned by the bare coupling 
constants where the geometry of spacctimc bears some resemblance with the 
kind of universe we observe around us. That is, the theory should create dy- 
namically an effective background geometry around which there are (small) 
quantum fluctuations. This is a very nontrivial property of the theory and one 
we are going to investigate in some detail. Computer simulations presented in 
these lectures confirm in a much more direct way the indirect evidence for such 
a scenario which we have known for some time and first reported in [11, 12]. 
They establish the de Sitter nature of the background spacetime, quantify the 
fluctuations around it, and set a physical scale for the universes we are dealing 
with. The main results of these investigations, without the numerical details, 
were announced in [13] and a detailed account of the results was presented in 
[14]. 

The remainder of these lecture notes is organized as follows: in Sec. 2 we 
describe the lattice formulation of four-dimensional quantum gravity. In Sec. 
3 the numerical results in 4d are summarized. We view these results as very 
important, but they also serve as a motivation for moving to two dimensions. 
While there is no propagating graviton in two-dimensional quantum grav- 
ity, it is a diffeomorphism- invariant theory and almost all of the conceptional 
problems mentioned above arc present there. Thus it is an important exercise 
to solve two-dimensional quantum gravity. Surprisingly this can be done in 
the lattice regularization known as "dynamical triangulation" . An important 
corollary is that (1) one can explicitly construct the continuum limit of the 
lattice theory and (2) show that it agrees with the so-called Liouville 2d quan- 
tum gravity theory. This latter theory is a continuum conformal field theory, 
explicit solvable, and, when viewed in the correct way, a diffeomorphism- 
invariant theory. Thus there is indeed no problem having a lattice regulariza- 
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tion of a diffeomorphism- invariant theory. In Sec. 4 we solve what is known as 
two-dimensional Euclidean quantum gravity. In Sec. 5 wc show how one can 
interpolate from Euclidean 2d quantum gravity to "Lorentzian" 2d quantum 
gravity which is a two-dimensional version of the four-dimensional gravity the- 
ory we have discussed in Sec. 3. Finally Sec. 7 discusses the results obtained 
and outlines perspectives. 

2 CDT 

The use of so-called causal dynamical triangulations (CDT) stands in the; 
tradition of [16], which advocated that in a gravitational path integral with 
the correct, Lorentzian signature of spacetime one should sum over causal 
geometries only. More specifically, we adopted this idea when it became clear 
that attempts to formulate a Euclidean non-perturbative quantum gravity 
theory run into trouble in spacetime dimension d larger than two as will be 
described below. 

This implies that we start from Lorentzian simplicial spacetimes with d = 4 
and insist that only causally well-behaved geometries appear in the (regular- 
ized) Lorentzian path integral. A crucial property of our explicit construction 
is that each of the configurations allows for a rotation to Euclidean signature. 
We rotate to a Euclidean regime in order to perform the sum over geometries 
(and rotate back again afterward if needed). We stress here that although 
the sum is performed over geometries with Euclidean signature, it is different 
from what one would obtain in a theory of quantum gravity based ab initio on 
Euclidean spacetimes. The reason is that not all Euclidean geometries with a 
given topology are included in the "causal" sum since in general they have no 
correspondence to a causal Lorentzian geometry. 

We refer to [6] for a detailed description of how to construct the class of 
piecewise linear geometries used in the Lorentzian path integral. The most 
important assumption is the existence of a global proper-time foliation. This 
is symbolically illustrated in Fig. 3 where we compare the construction to the 
one of ordinary quantum mechanics: the path integral of ordinary quantum 
mechanics is regularized as a sum over piecewise linear paths from point a;, 
to point Xf in time tf — U. The time steps have length a and the continuum 
limit is obtained when the length a of these "building blocks" goes to zero. 
Similarly, in the quantum gravity case we have a sum over four-geometries, 
"stretching" between two three-geometries separated a proper time t and con- 
structed from four-dimensional building blocks, as described below. On the 
figure we show for illustrational simplicity only a single space-time history 
and replace the three-dimensional spatial geometries with a one-dimensional 
one of S'^-topology. Moving from left to right we have a time foliation where 
at each discrete time step space is represented by a circle. Neighbouring cir- 
cles are then connected by piecewise flat biiilding blod^s, usually triangles, as 
illustrated in Fig. 17 in Sec. 5. In the "real" four-dimensional case, the spatial 




Fig. 3. Piecewise linear space-time liistories in quantum mechanics and in (1+1)- 
dimensional quantum gravity. In tlie gravity case we siiow only a single space-time 
history, while in the quantum particle case we show many such histories as well as 
the average path (thick line). 



slices of topology will be replaced by spatial slices of topology S^, and 
neighbouring spatial slices are then connected by four-simplices as illus- 
trated in Fig. 4 and described in detail below. There is an important difference 
between the quantum mechanical sum over paths and our sum over geome- 
tries with a time foliation: the time t in the quantum mechanical example is an 
external parameter, while the time t in the case of quantum gravity is intrin- 
sic. Also, again for the purpose of illustration, the two-dimensional geometry 
has been drawn embedded in three-dimensional space, but in the path inte- 
gral implementing the summation over geometries there is no such embedding 
present. Finally, we cannot refrain from mentioning that the paths shown for 
the quantum mechanical particle are in fact typical paths which appear for 
the (Euclideanized) path integral of a particle placed in an external potential. 
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They are picked out from an actual Monte Carlo simulation of such a physical 

system. Similarly, the 2d surface is a surface picked out from a Monte Carlo 
simulation of 2d quantum gravity and thus corresponds to a typical 2d surface 
which appears in the path integral. This is the reason for the somewhat poor 
graphic representation: there is no natural length-preserving representation 
of the surface in 3d such that the surface is not self-intersecting. For better 
graphic illustrations (animations) of the 2d surfaces which appear in the 2d 
quantum-gravitational path integral we refer to the link [7]. 

As mentioned above, we assume that the spacetime topology is that of 
X R, the spatial topology being that of merely for convenience. The 
spatial geometry at each discrete proper-time step f„ is represented by a 
triangulation of S^, made up of equilateral spatial tetrahedra with squared 
side-length = > 0. In general, the number Nz{tn) of tetrahedra and how 
they are glued together to form a picccwise flat three-dimensional manifold 
will vary with each time-step tn- In order to obtain a four-dimensional tri- 
angulation, the individual three-dimensional slices must still be connected in 
a causal way, preserving the S'^-topology at all intermediate times t between 
tn and tn+i ^- This is done as illustrated in Fig. 4, introducing what we call 
(4,l)-simplices and (3,2)-simplices. More precisely, a (4, l)-simplex is a four- 
simplex with foiu of its vertices (i.e. a boundary tetrahedron) belonging to the 
triangulation of S^{tn), the time-slice corresponding to time tn, and the fifth 
vertex belonging to the triangulation of S^{tn+i), the time-slice correspond- 
ing to time Similarly, a (3,2) simplex has three vertices, i.e. a triangle, 
belonging to the triangulation of S^{tn) and two vertices, i.e. a link, belong- 
ing to the triangulation of S^{tn+i)- We have also simplices of type (1,4) and 
(2,3), which are defined in an obvious way, interchanging the role of S^{tn) 
and S^{tn+i)- One can show that two triangulations of S^{tn) and S^{tn+i) 
can be "connected" by these four building blocks glued together in a suitable 
way such that we have a four-dimensional triangulation of x [0,1]. Also, 
two given triangulations of S^{tn) and S^{tn+i) can be connected in many 
ways compatible with the topology 5^ x [0, 1]. In the path integral we will be 
summing over all possible ways to connect a given triangulation S^{tn) to a 
given triangulation of S^{tn+i) compatible with the topology x [0,1]. In 
addition we will sum over all 3d triangulations of 5^ at all times tn- 

We allow for an asymmetry between temporal and spatial lattice length 
assignments. Denote by it and £g the length of the time-like links and the 
space- like links, respectively. Then f| = —ai"^, a > 0. The explicit rotation 
to Euclidean signature is done by performing the rotation a —a in the 
complex lower half-plane, |a| > 7/12, such that we have = \a\£l (see [6] for 
a discussion). 



^ This implies the absence of branching of the spatial universe into several dis- 
connected pieces, so-called baby universes, which (in Lorcntzian signature) would 
inevitably be associated with causality violations in the form of degeneracies in 
the light cone structure, as has been discussed elsewhere (see, for example, [15]). 
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(4,1) (3,2) 

Fig. 4. (4,1) and a (3,2) simplices connecting two neighbouring spatial slices. We 
also have symmetric (1,4) and (2,3) simplices with a vertex and a line, respectively, 
at time t and a tetrahedron and a triangle, respectively, at time t+1. For simplicity 
we denote the total number of (4,1) and (1,4) simplices by A'^'''^' and similarly the 
total number of (3,2) and (2,3) simplices by N^'^ . 

The Einstein-Hilbert action has a natural geometric implementation 
on piecewise linear geometries in the form of the Regge action. This is given 
by the sum of the so-called deficit angles around the two-dimensional "hinges" 
(subsimplices in the form of triangles) , each multiplied with the volume of the 
corresponding hinge. In view of the fact that we are dealing with piecewise 
linear, and not smooth metrics, there is no unique "approximation" to the 
usual Einstein-Hilbert action, and one could in principle work with a different 
form of the gravitational action. We will stick with the Regge action, which 
takes on a very simple form in our case, where the piecewise linear manifold 
is constructed from just two different types of building blocks. After rotation 
to Euclidean signature one obtains for the action (see [8] for details) 

sf = / '^'^Vai-R + ^ (12) 

Sfss'' = + 6A)No + «;4(^f + A^f ■'^) + zi(2A^f + Tvf '")), 

where Nq denotes the total number of vertices in the four-dimensional tri- 
angulation and N^'^'^^ and n[^'^^ denote the total number of the four- 
simplices described above, so that the total number N4 of four-simplices is 
N4 = N^'^'^'^ + n[^'^\ The dimensionless coupling constants kq and K4 are 
related to the bare gravitational and bare cosmological coupling constants, 
with appropriate powers of the lattice spacing a already absorbed into kq a-nd 
K4. The asymmetry parameter A is related to the parameter a introduced 
above, which describes the relative scale between the (squared) lengths of 
space- and time-like links. It is both convenient and natural to keep track of 
this parameter in our set-up, which from the outset is not isotropic in time 
and space directions, see again [8] for a detailed discussion. Since we will in 
the following work with the path integral after Wick rotation, let us redefine 
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a := —a [8], which is positive in the Euclidean domain.^ For future reference, 
the Euclidean four-volume of our universe for a given choice of a is given by 

= (74(0 «^ ^f'^ = <54(0 ^4/(1 + 0, (13) 

where ^ is the ratio 

^ = 7Vf2)/Arfi), (14) 

and Ci{^) a* is a measure of the "effective four- volume" of an "average" four- 
simplex. ^ will depend on the choice of coupling constants in a rather compli- 
cated way (for a detailed discussion we refer to [6, 14]). 

The path integral or partition function for the CDT version of quantum 
gravity is now 

Z(G,yl) = |p[ff]e-^i"[^] - Z(«o,«4,4) = E^<^"'"^^'' (15) 

where the summation is over all causal triangulations T of the kind described 
above, and we have dropped the superscript "Regge" on the discretized action. 
The factor l/CV is a symmetry factor, given by the order of the automor- 
phism group of the triangulation T. The actual set-up for the simulations is 
as follows. We choose a fixed number of spatial slices at proper times t\, 
t2 = ti + at, up to In = ti + {N— l)at, where At = at is the discrete lattice 
spacing in temporal direction and T = Nat the total extension of the universe 
in proper time. For convenience we identify %+i with ti, in this way imposing 
the topology x 5^ rather than I x S^. This choice does not affect physical 
results, as will become clear in due course. 

Our next task is to evaluate the non-perturbative sum in (15), if possible, 
analytically. This can be done in spacetime dimension d = 2 ([10, 17] (and 
we discuss this in detail below) and at least partially in d = 3 [18, 19], but 
presently an analytic solution in four dimensions is out of reach. However, we 
are in the fortunate situation that Z{ko, m, A) can be studied quantitatively 
with the help of Monte Carlo simulations. The type of algorithm needed to up- 
date the piecewise linear geometries has been around for a while, starting from 
the use of dynamical triangulations in bosonic string theory (two-dimensional 
Euclidean triangulations) [20, 21, 22] and was later extended to their ap- 
plication in Euclidean four-dimensional quantum gravity [23, 24]. In [6] the 
algorithm was modified to accommodate the geometries of the CDT set-up. 
The algorithm is such that it takes the symmetry factor Cr into account 
automatically. 

We have performed extensive Monte Carlo simulations of the partition 
function Z for a number of values of the bare coupling constants. As reported 
in [8] , there are regions of the coupling constant space which do not appear 



^ The most symmetric choice is a = 1, corresponding to vanishing asymmetry, 
A = Q. 
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Fig. 5. The phases A, B and C in the coupling constant space (kq, ^, K4)- Phase C 
is the one where extended four-dimensional geometries emerge. 

relevant for continuum physics in that they seem to suffer from problems sim- 
ilar to the ones found earlier in Euclidean quantum gravity constructed in 
terms of dynamical triangulations, which essentially led to its abandonment 
in d > 2. What is observed in Euclidean four-dimensional quantum gravity is 
the following: when the (inverse, bare) gravitational coupling kq is sufficiently 
large one sees so-called branched polymers, i.e. not really a four-dimensional 
universe, but a universe which branches out like a tree with so many branches 
that it becomes truly fractal when the number of four-simplices becomes infi- 
nite, and its Hausdorff dimension is 2. Such triangulations represent the most 
extended triangulations one can construct unless one explicitly forbids branch- 
ing. When the (inverse, bare) gravitational coupling kq is sufficiently small one 
observes a totally crumpled universe with almost no extension. In this phase 
there exist vertices of very high order and the connectivity of the triangulation 
is such that it is possible to move from any four-simplex to any other crossing 
only a few neighbouring four-simplices. The Hausdorff dimension of such a 
triangulation is infinite in the limit where the number of four-simplices goes 
to infinity. These two phases, the crumpled and the branched polymer phase, 
are separated by a phase transition line along which there is a first-order tran- 
sition. It was originally hoped that one could find a point on the critical line 
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where the first-order transition becomes second order and which could then 

be used as a fixed point where a continuum theory of quantum gravity could 
be defined along the lines suggested by eqs. (10) and (11). However, such a 
second-order transition was not found [9], and eventuallythe idea of a theory of 
four-dimensional Euclidean quantum gravity was abandoned. A new principle 
for selecting the class of geometries one should use in the path integral was 
needed and this led to the suggestion to include only causal triangulations in 
the sum over spacctimc histories. 

When we include only the causal triangulations in the path integral, we 
still see a remnant of the Euclidean structure just described, namely, when 
the (inverse, bare) gravitational coupling kq is sufficiently large, the Monte 
Carlo simulations exhibit a sequence in time direction of small, disconnected 
universes, none of them showing any sign of the scaling one would expect 
from a macroscopic universe. We denote this phase by A. We believe that this 
phase of the system is a Lorentzian version of the branched-polymer phase 
of Euclidean quantum gravity. By contrast, when A is sufficiently small, the 
simulations reveal a universe with a vanishing temporal extension of only 
a few lattice spacings, ending both in past and future in a vertex of very 
high order, connected to a large fraction of all vertices. This phase is most 
likely related to the so-called crumpled phase of Euclidean quantum gravity. 
We denote this phase by B. The crucial and new feature of the quantum 
superposition in terms of causal dynamical triangulations is the appearance of 
a region in coupling constant space which is different and interesting and where 
continuum physics may emerge. It is in this region that we have performed the 
simulations discussed here and where work up to now has already uncovered 
a number of intriguing physical results [11, 12, 8, 25]. In Fig. 6 we have shown 
how different configurations look in the three phases discussed above, and in 
Fig. 5 we have shown the tentative phase diagram in the coupling constant 
space of Ko, K4 and A. A "critical" surface is shown in the figure. Keeping kq 
and A fixed, K4 acts as a chemical potential for N4; the smaller K4, the larger 
(A^4). At some critical value K4{ko, A), depending on the choice of kq and A, 
{N4) 00. For K4 < K4{k(),A) the partition function is plainly divergent 
and not defined. When we talk about phase transitions we are always at the 
"critical" surface 

K4 = K4{kq,A), (16) 

simply because we cannot have a phase transition unless A'4 = 00. We put 
"critical" into quotation marks since it only means that we probe infinite 
four- volume. No continuum limit is necessarily associated with a point on this 
surface. To decide this issue requires additional investigation. A good anal- 
ogy is the Ising model on a finite lattice. To have a genuine phase transition 
for the Ising model we have to take the lattice volume to infinity since there 
arc no genuine phase transitions for finite systems. However, just taking the 
lattice volume to infinity is not sufficient to ensure critical behaviour of the 
Ising model. We also have to tune the coupling constant to its critical value. 
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Fig. 6. Typical configurations in the phases A, B and C (lowest figure). Phase C is 
the one where extended four-dimensional geometries emerge. 



Being on the "critical" surface, or rather "infinite- volume" surface (16), we 
can discuss various phases, and these are the ones indicated in the figure. The 
different phases are separated by phase transitions, vifhich might be first-order. 
However, we have not yet conducted a systematic investigation of the order 
of the transitions. Looking at Fig. 5, we have two hues of phase transitions, 
separating phase A and phase C and separating phase B and phase C re- 
spectively. They meet in the point indicated on the figure. It is tempting to 
speculate that this point might be associated with a higher-order transition, 
as is common for statistical systems in such a situation. We will return to this 
point later. 

In the Euclideanized setting the value of the cosmological constant deter- 
mines the spacetime volume V4 since the two appear in the action as conjugate 
variables. We therefore have (V4) ~ Gj A in a continuum notation, where G is 
the gravitational coupling constant and A the cosmological constant. In the 
computer simulations it is more convenient to keep the four-volume fixed or 
partially fixed. We will implement this by fixing the total number of four- 
simplices of type Nf"'^^ or, equivalently, the total number of tetrahedra 
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making up the spatial triangulations at times ti, i = 1, . . . ,N, 

N, = f2N,{t,) = ^Ni'''\ (17) 

i=l 

We know from the simulations that in the phase of interest {N^'^'^'') oc (A'^"^'^-') 
as the total volume is varied [8]. This effectively implies that we only have 
two bare coupling constants kq, ^ in (15), while we compensate by hand for 
the coupling constant K4 by studying the partition function Z{ko, A; N^'^'^^) 
for various Nji'^'^\ To keep track of the ratio ^{kq, A) between the expectation 
value {N^^'^^} and N^^'^\ which depends weakly on the coupling constants, 
we write (c.f. eq. (14)) 

(iV4) = ivf + (ivf ^) = A^f ')(1 + e(«o, A)). (18) 

For all practical purposes we can regard A'4 in a Monte Carlo simulation as 
fixed. The relation between the partition function we use and the partition 
function with variable four- volume is given by the Laplace transformation 

Z{ko, K4, A)= I dNi e-'"'^' Z{kq, N4, A), (19) 
Jo 

where strictly speaking the integration over N4 should be replaced by a sum- 
mation over the discrete values A^4 can take. Returning to Fig. 5, keeping A'4 
fixed rather than fine-tuning K4 to the critical value k'I implies that one is 
already on the "critical" surface drawn in Fig. 5, assuming tliat A'4 is siiffi- 
cicntly large (in principle infinite). Whether N4 is sufficiently large to quafify 
as "infinite" can be investigated by performing the computer simulations for 
different 7V4's and comparing the results. This is a technique we will use over 
and over again in the following. 

3 Numerical results 

The Monte Carlo simulations referred to above will generate a sequence of 

spacetime histories. An individual spacetime history is not an observable, in 
the same way as a path x{t) of a particle in the quantum-mechanical path 
integral is not. However, it is perfectly legitimate to talk about the expectation 
value {x{t)) as well as the fluctuations around {x{t)). Both of these quanti- 
ties are in principle calculable in quantum mechanics. Let us make a slight 
digression and discuss this in some detail since it illustrates well the picture 
we also hope emerges in a theory of quantum gravity. Consider the particle 
example shown in Fig. 3. We have a particle moving from Xi at ti to Xf at 
tf. In general there will be a classical motion of the particle satisfying these 
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boundary conditions (we will assume that for simplicity). If H can be consid- 
ered small compared to the other parameters entering into the description of 
the system, the classical path will be a good approximation to {x{t)) accord- 
ing to Ehrenfest's theorem. In Fig. 3 the smooth curve represents {x{t)). In 
the path integral we smn over all continuous paths from {xiAi) to (xf,tf) as 
illustrated in Fig. 3. However, when all other parameters in the problem are 
large compared to H we expect a "typical" path to be close to {x{t)) which 
also will be close to the classical path. Let us make this explicit in the simple 
case of the harmonic oscillator. Let Xd {t) denote the solution to the classical 
equations of motion such that XdiU) = Xi and Xci{tf) = Xf. For the harmonic 
oscillator the decomposition 



leads to an exact factorization of the path integral thanks to the quadratic 
nature of the action. The part involving Xci{t) gives precisely the classical 
action and the part involving y(t) the contributions from the fluctuations, 
independent of the classical part. Taking the classical path to be macroscopic 
gives a picture of a macroscopic path dressed with small quantum fluctuations, 
small because they are independent of the classical motion. Explicitly we have 
for the fluctuations (Euclidean calculation): 



Thus the harmonic oscillator is a simple example of what we hope for in 

quantmn gravity: Let the size of the system be macroscopic, i.e. Xci{t) is 
macroscopic (put in by hand): then the quantum fluctuations around this 
path are small and of the order 



We hope this translates into the description of our universe: the macroscopic 
size of the universe dictated by the (inverse) cosmological constant in any 
Euclidean description (trivial to show in the model by simply differentiating 
the partition function with respect to the cosmological constant and in the 
simulations thus put in by hand) and the small quantum fluctuations dictated 
by the other coupling constant, namely, the gravitational coupling constant. 

3.1 The emergent de Sitter background 

Obviously, there are many more dynamical variables in quantum gravity than 
there are in the particle case. We can still imitate the quantum-mechanical 

situation by picking out a particular one, for example, the spatial three- volume 
V3{t) at proper time t. We can measure both its expectation value (V3(t)) as 



x{t) = Xci{t) + y{t), y{ti) = y{tf) = 
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well as fluctuations around it. The former gives us information about the 
large-scale "shape" of the universe we have created in the computer. First we 
will describe the measurements of (V3(t)), keeping a more detailed discussion 
of the fluctuations to Sec. 3.2 below. 

A "mcasmement" of V3(t) consists of a table iV3(i), where i — 1, ... ,7V 
denotes the number of time-slices. Recall from Sec. 2 that the sum over slices 
'Yli'i=i ^3{i) is kept constant. The time axis has a total length of N time steps, 
where TV = 80 in the actual simulations, and we have cyclically identifled 
time-slice N + 1 with time-slice 1. 

What we observe in the simulations is that for the range of discrete vol- 
umes N4 under study the universe does not extend (i.e. has appreciable three- 
volume) over the entire time axis, but rather is localized in a region much 
shorter than 80 time slices. Outside this region the spatial extension Ns{i) 
will be minimal, consisting of the minimal number (five) of tetrahedra needed 
to form a three-sphere S^, plus occasionally a few more tetrahedra."^ This thin 
"stalk" therefore carries little four-volume and in a given simulation we can 
for most practical purposes consider the total four-volume of the remainder, 
the extended universe, as fixed. 

In order to perform a meaningful average over geometries which explicitly 
refers to the extended part of the universe, we have to remove the transla- 
tional zero mode which is present. We refer to [14] for a discussion of the 
procedure. Having defined the centre of volume along the time-direction of 
our spaeetimc configurations we can now perform superpositions of configura- 
tions and define the average (iV3(«)) as a function of the discrete time i. The 
results of measuring the average discrete spatial size of the universe at various 
discrete times i are illustrated in Fig. 7 and can be succinctly summarized by 
the formula 

where N3{i) denotes the number of three-simplices in the spatial slice at dis- 
cretized time i and A^4 the total number of four-simplices in the entire universe. 
Since we are keeping A^j^'^-* fixed in the simulations and since ^ changes with 
the choice of bare coupling constants, it is sometimes convenient to rewrite 
(20) as 

m\i) = lAff 1) ^ i- cos3 f ^- ^ , (21) 

where §0 is defined by so(^4^'^^)^^^ = soN\^^. Of course, formula (20) is only 
valid in the extended part of the universe where the spatial three-volumes are 

larger than the minimal cut-ofF size. 

This kinematic constraint ensures that the triangulation remains a simplicial 
manifold in which, for example, two d-simplices are not allowed to have more 
than one (d — l)-simplex in common. 
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Fig. 7. Background geometry {N'i{i)): MC measurements for fixed A^'^'''^-' — 160.000 
{Ni — 362.000) and best fit (20) yield indistinguishable curves at given plot resolu- 
tion. The bars indicate the average size of quantum fluctuations. 

The data shown in Fig. 7 have been collected at the particular values 
{ko,A) = (2.2,0.6) of the bare coupling constants and for N4 — 362.000 
(corresponding to ivj'*'^'' = 160.000). For this value of (kq, ^) we have verified 
relation (20) for N4 ranging from 45.500 to 362.000 building blocks (45.500, 
91.000, 181.000 and 362.000). After rescaling the time and volume variables 
by suitable powers of A^4 according to relation (20) , and plotting them in the 
same way as in Fig. 7, one finds almost total agreement between the curves for 
different spacetime volumes. This is illustrated in Fig. 8. Thus we have here 
a beautiful example of finite-size scaling, and at least when we discuss the 
average three- volume V^^t) all our discretized volumes are large enough 
that we can treat them as infinite, in the sense that no further change will 
occur for larger 7V4. 

By contrast, the quantum fluctuations indicated in Fig. 7 as vertical bars 
are volume-dependent and will be the larger the smaller the total four-volume, 
see Sec. 3.2 below for details, eq. (20) shows that spatial volumes scale accord- 
ing to N^^'^ and time intervals according to N^^'^, as one would expect for a 
genuinely /owr-dimensional spacetime and this is exactly the scaling we have 
used in Fig. 8. This strongly suggests a translation of (20) to a continuum 
notation. The most natural identification is given by 

VmV,''(t) = v,^cos'(^^y (22) 

where we have made the identifications 
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I = ^171, ^U^t VsiU) = 2CiN3{i)a\ (23) 
such that we have 

' dt,/grtV3{t) = Vi. (24) 



In (23), is the constant proportionality factor between the time t and 

genuine contimmm proper time t, r = ^/gtl t- (The combination Ati^fg^V^ 
contains Ci, related to the four- volume of a four-simplex rather than the three- 
volume corresponding to a tetrahedron, because its time integral must equal 
Vi). Writing V4 = 87r^i?'*/3, and ,Jgtt = R/B, cq. (22) is seen to describe a 
Euclidean de Sitter universe (a four-sphere, the maximally symmetric space 
for positive cosmological constant) as our searched-for, dynamically gener- 
ated background geometry! In the parametrization of (22) this is the classical 
solution to the action 

whore fc2 = 9(27r^)^/'' and A is a Lagrange multiplier, fixed by requiring that 
the total four- volume be V4, J dt^/gu Vzit) = V4. Up to an overall sign, this is 
precisely the Einstein-Hilbert action for the scale factor a{t) of a homogeneous, 
isotropic universe (rewritten in terms of the spatial three- volume ¥3(1) = 
2TT^a{t)^), although we of course never put any such simplifying symmetry 
assumptions into the CDT model. 

A discretized, dimensionless version of (25) is 



^discr 



^1 2^ [ J^) + ^^ATg' {i) ] , (26) 



where /c2 oc ^2- This can be seen by applying the scaling (20), namely, A^3(i) = 
-^4^*'^3(si) and Si = i/N^'^. This enables us to finally conclude that the 
identifications (23) when used in the action (26) lead naively to the continuum 
expression (25) under the identification 



^2 , /r< s2 

^^aVC4So_ (27) 
fci 3v^ 

Next, let us comment on the universality of these results. First, we have 
checked that they are not dependent on the particular definition of time- 
slicing we have been using, in the following sense. By construction of the 
piecewise linear CDT-geometries we have at each integer time step ti = iat a. 
spatial surface consisting of A^3(i) tetrahedra. Alternatively, one can choose as 
reference slices for the measurements of the spatial vohimo non-integer values 
of time, for example, all time slices at discrete times i — 1/2, i = 1,2,... . In 
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Fig. 8. Rescaling of time and volume variables according to relation (20) for A'4 = 
45.500, 91.000, 181.000 and 362.000. The plot also include the curve (20). More 
precisely: a oc i/N^^"^ and P{cj) oc N:i{i)/N^^^ . 

this case the "triangulation" of the spatial three-spheres consists of tetrahedra 
- from cutting a (4,1)- or a (l,4)-simplex half-way - and "boxes", obtained 
by cutting a (2,3)- or (3,2)-simplex (the geometry of this is worked out in 
[26] ) . We again find a relation hke (20) if we use the total number of spatial 
building blocks in the intermediate slices (tetrahedra-|-boxes) instead of just 
the tetrahedra. 

Second, we have repeated the measurements for other values of the bare 
coupling constants. As long as we stay in the phase where an extended uni- 
verse is observed, the phase C in Fig. 5, a relation like (20) remains valid. In 
addition, the value of Sq, defined in eq. (20), is almost unchanged until we 
get close to the phase transition lines beyond which the extended universe 
disappears. Only for the values of kq around 3.6 and larger will the measured 
{N3{t)) differ significantly from the value at 2.2. For values larger than 3.8 
(at A — 0.6), the universe will disintegrate into a number of small and dis- 
connected components distributed randomly along the time axis, and one can 
no longer fit the distribution {N3{t)) to the formula (20). Later we will show 
that while sq is almost unchanged, the constant ki in (26), which governs 
the quantum fluctuations around the mean value {N3{t)), is more sensitive to 
a change of the bare coupling constants, in particular, in the case where we 
change kq (while leaving A fixed). 
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3.2 Fluctuations around de Sitter space 

In the following we will test in more detail how well the actions (25) and 
(26) describe the computer data. A crucial test is how well it describes the 

quantum fluctuations around the emergent de Sitter background. 
The correlation function (the covariance matrix C) is defined by 

CNAi,i') = {6N3{t)SN3{i')), SNaii) = Naii) - N^ii), (28) 

where we have included an additional subscript A^4 to emphasize that N4 is 
kept constant in a given simulation. 

The first observation extracted from the Monte Carlo simulations is that 
under a change in the four- volume Cjv4 («,«') scales as'' 

Cn, {i, i') = NiF {ilNl'\ i'/Nl'^) , (29) 

where F is a universal scaling function. This is illustrated by Fig. 9 for the 
rescaled version of the diagonal part C^^(i, i), corresponding precisely to the 
quantum fluctuations {{5N3{i))'^)^/'^ of Fig. 7. While the height of the curve in 
Fig. 7 will grow as ^^4^^, the superimposed fluctuations will only grow as N^^"^ . 
We conclude that for fixed hare coupling constants the relative fluctuations will 
go to zero in the infinite-volume limit . 

Let us rewrite the minisuperspace action (25) for a fixed, finite four-volume 

Vi in terms of dimensionless variables by introducing s = t/V^''^ and Vait) = 
V!'\,{s): 

now assuming that / ds^/g^ ^^{s) — 1, and with g^s = gtt - The same rewriting 
can be done to (26) which becomes 

Saiscr = fci^E^^ [-^^ j^ n3(..+0-n3(..) y_^-^^v3^^^^j ^ ^3^^ 

where A^3(i) = A^4^^n3(sj) and Si = i/N^^'^. 

Prom the way the factor \/Ni appears as an overall scale in eq. (31) it is 
clear that to the extent a quadratic expansion around the effective background 
geometry is valid one will have a scaling 

{5N^{i)6Nr,{i')) = Nl'^{5m{ti)5m{ti,)) = NiF{ti,ti,), (32) 

^ We stress again that the form (29) is only valid in that part of the universe whose 
spatial extension is considerably larger than the minimal 5"^ constructed from 5 
tetrahedra. (The spatial volume of the stalk typically fluctuates between 5 and 
15 tetrahedra.) 




where ti = i/N^^^. This imphes that (29) provides additional evidence for the 
validity of the quadratic approximation and the fact that our choice of action 
(26), with ki independent of is indeed consistent. 

To demonstrate in detail that the fuU function F{t,t') and not only its 
diagonal part is described by the effective actions (25), (26), let us for conve- 
nience adopt a continuum language and compute its expected behaviour. Ex- 
panding (25) around the classical solution according to V3{t) = V^'-{t) + x{t), 
the quadratic fluctuations are given by 

w.wn)=/i>.(.).(w*>-u/-'.<.)-<...'W) 

= */-■((,('), (33) 

where 'Dx{s) is the normalized measure and the quadratic form M{t,t') is 
determined by expanding the effective action S to second order in x{t)^ 

S{Vs) = S{V,^') + J dt x{t)Hx{t)- (34) 



In expression (34), H denotes the Hermitian operator 
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d 1 d 4 

dicos3(i/-B)dt B^cos^{t/B)' ^ ' 

which must be diagonalized under the constraint that J dt,/gu x{t) = 0, since 
V4 is kept constant. 

Let e^"^ (i) be the eigenfunctions of the quadratic form given by (34) with 
the volume constraint enforced, ordered according to increasing eigenvalues 
A„. As we will discuss shortly, the lowest eigenvalue is Ai = 0, associated 
with translational invariance in time direction, and should be left out when 
we invert M{t, t'), because we precisely fix the centre of volume when making 
our measurements. Its dynamics is therefore not accounted for in the correlator 
C{t,t'). 

If this cosmological continuum model were to give the correct description 
of the computer-generated universe, the matrix 
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e(i)r 



M-\t,t') = y2 4 — (36) 

should be proportional to the measured correlator C{t,t'). Fig. 10 shows the 
eigcnfunctions e'^^^(t) and e'^^^(t) (with two and four zeros respectively), calcu- 
lated from H with the constraint / dt^Jgu x{t) = imposed. Simultaneously 
we show the corresponding eigenfunctions calculated from the data, i.e. from 
the matrix C{t,t'), which correspond to the (normalizablc) eigcnfunctions 
with the highest and third-highest eigenvalues. The agreement is very good, 
in particular, when taking into consideration that no parameter has been ad- 
justed in the action (we simply take B = soNl^^At in (22) and (34), which 
gives B = UA7at for Ni = 362.000). 

The reader may wonder why the first eigenfunction exhibited has two zeros. 
As one would expect, the ground state eigenfunction e^^\t) of the Hamiltonian 
(35), corresponding to the lowest eigenvalue, has no zeros, but it does not 
satisfy the volume constraint / dty/gu x{t) = 0. The eigenfunction e^^\t) of 
H with next-lowest eigenvalue has one zero and is given by the simple analytic 
function 

where c is a constant. One realizes immediately that e^"'^'' is the translational 
zero mode of the classical solution V^\t) (oc cos^ t/B). Since the action is 
invariant under time translations, we have 

S{V,^\t + At)) = SiV,^\t)), (38) 

and since V^^t) is a solution to the classical equations of motion we find to 
second order (using the definition (37)) 

S{V,^'{t + At)) = S{V,^\t)) + / dt e^'\t)M^\t), (39) 

consistent with e(^^(t) having eigenvalue zero. 

It is clear from Fig. 10 that some of the eigenfunctions of H (with the 
volume constraint imposed) agree very well with the measured eigenfunctions. 
All even eigenfunctions (those symmetric with respect to reflection about 
the symmetry axis located at the centre of volume) turn out to agree very 
well. The odd eigenfunctions of H agree less well with the eigenfunctions 
calculated from the measured C(t,t'). The reason seems to be that we have 
not managed to eliminate the motion of the centre of volume completely from 
our measurements. There is an inherent ambiguity in fixing the centre of 
volume of one lattice spacing, which turns out to be sufficient to reintroduce 
the zero mode in the data. Suppose we had by mistake misplaced the centre 
of volume by a small distance At. This would introduce a modification 

AV. = '-Ym At (40) 
dt ^ ' 
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proportional to the zero mode of the potential V^\t). It follows that the zero 
mode can re-enter whenever we have an ambiguity in the position of the centre 
of volume. In fact, we have found that the first odd eigenfunction extracted 
from the data can be perfectly described by a linear combination of e'^^ {t) and 
e^^\t). It may be surprising at first that an ambiguity of one lattice spacing 
can introduce a significant mixing. However, if we translate AV^ from eq. (40) 
to "discretized" dimensionless units using ^3(1) ~ N^^^ cosii/Nl^"^), we find 
that AV3 ~ which because of {{SN:}{i))^) ~ is of the same order of 

magnitude as the fluctuations themselves. In our case, this apparently does 
affect the odd eigenfunctions. 

One can also compare the data and the matrix M~^(t,t') calculated from 
(36) directly. This is illustrated in Fig. 11, where we have restricted ourselves 
to data from inside the extended part of the universe. We imitate the con- 
struction (36) for M~^, using the data to calculate the eigenfunctions, rather 
than H. One could also have used C{t,t') directly, but the use of the eigen- 
functions makes it somewhat easier to perform the restriction to the bulk. 
The agreement is again good (better than 15% at any point on the plot), 
although less spectacular than in Fig. 10 because of the contribution of the 
odd eigenfunctions to the data. 



3.3 The size of the universe and the flow of G 

It is natural to view the coupling constant G in front of the effective action 
for the scale factor as the gravitational coupling constant G. The effective 
action which described our computer-generated data was given by eq. (25) 
and its dimensionless lattice version by (26). The computer data allows us to 
extract k\ oc a^/G, a being the spatial lattice spacing, the precise constant of 
proportionality being given by eq. (27): 




(41) 

For the bare coupling constants {ko,A) = (2.2,0.6) we have high-statistics 

measurements for A^4 ranging from 45.500 to 362.000 four-simplices (equiva- 

lently, A^j^'^'' ranging from 20.000 to 160.000 four-simplices). The choice of A 
determines the asymmetry parameter a, and the choice of {kq, A) determines 
the ratio ^ between TVp'^-* and N^'^'^\ This in turn determines the "effective" 
four-volume G4 of an average four-simplex, which also appears in (41). The 
number §0 in (41) is determined directly from the time extension Tuniv of the 
extended universe according to 

Tuniv = 7r so (TVf'))'^'. (42) 

Finally, from our measurements we have determined ki = 0.038. Taking ev- 
erything together according to (41), we obtain G w 0.23a^ or £pi w 0.48a, 
where ipi = VG is the Planck length. 
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Fig. 11. Comparing data for the extended part of the universe: measured C{t,t') 
(above) versus M~^{t,t') obtained from analytical calculation (below). The agree- 
ment is good, and would have been even better had we included only the even 
modes. 

From the identification of the volume of the four-sphere, V4 = 87r^i?*/3 — 
C^N^^'^"^ a^, we obtain that R — 3.1a. In other words, the linear size irR 
of the quantum de Sitter universes studied here lies in the range of 12-21 
Planck lengths for A'4 in the range mentioned above and for the hare coupling 
constants chosen as {kq^A) — (2.2,0.6). 

Our dynamically generated universes are therefore not very big, and the 
quantum fluctuations around their average shape are large as is apparent from 
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Fig. 7. It is rather surprising that the semi-classical minisuperspace formu- 
lation is applicable for universes of such a small size, a fact that should be 
welcome news to anyone performing semi-classical calculations to describe the 
behaviour of the early universe. However, in a certain sense our lattices are 
still coarse compared to the Planck scale £pi because the Planck length is 
roughly half a lattice spacing. If we are after a theory of quantum gravity 
valid on all scales, we are in particular interested in uncovering phenomena 
associated with Planck-scalc physics. In order to collect data free from un- 
physical short-distance lattice artifacts at this scale, we would ideally like to 
work with a lattice spacing much smaller than the Planck length, while still 
being able to set by hand the physical volume of the universe studied on the 
computer. 

The way to achieve this, under the assumption that the coupling constant 
G of formula (41) is indeed a true measure of the gravitational coupling con- 
stant, is as follows. We are free to vary the discrete four- volume A^4 and the 
bare coupling constants (acq, ^) of the Regge action (see [8] for further details 
on the latter). Assuming for the moment that the semi-classical minisuper- 
space action is valid, the effective coupling constant ki in front of it will be 
a function of the bare coupling constants {kq, A), and can in principle be de- 
termined as described above for the case (kq.A) = (2.2,0.6). If we adjusted 
the bare coupling constants such that in the limit as N4 00 both 

Vi^A^V and G /ki{Ko, A) (43) 

remained constant (i.e. fci(Ko, A) ~ l/y/Ni), we would eventually reach a re- 
gion where the Planck length was significantly smaller than the lattice spacing 
a, in which event the lattice could be used to approximate spacetime structures 
of Planckian size and we could initiate a genuine study of the sub-Planckian 
regime. Since we have no control over the effective coupling constant fci, the 
first obvious question which arises is whether we can at all adjust the bare cou- 
pling constants in such a way that at large scales we still see a four-dimensional 
universe, with ki going to zero at the same time. The answer seems to be in 
the affirmative, as we will go on to explain. 

Fig. 12 shows the results of extracting fci for a range of bare coupling 
constants for which we still observe an extended universe. In the top figure A = 
0.6 is kept constant while kq is varied. For kq sufficiently large we eventually 
reach a point where a phase transition takes place (the point in the square 
in the bottom right-hand corner is the measurement closest to the transition 
we have looked at). For even larger values of kq, beyond this transition, the 
universe disintegrates into a number of small universes, in a CDT-analogue of 
the branched-polymer phase of Euclidean quantum gravity. The plot shows 
that the effective coupling constant ki becomes smaller and possibly goes to 
zero as the phase transition point is approached, although our current data 
do not yet allow us to conclude that ki does indeed vanish at the transition 
point. 
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Fig. 12. The measured effective coupling constant ki as function of the bare kq (top, 
A — 0.6 fixed) and the asymmetry A (bottom, kq = 2.2 fixed). The marked point 
near the middle of the data points sampled is the point {no, A) — (2.2,0.6) where 
most measurements in the remainder of the paper were taken. The other marked 
points are those closest to the two phase transitions, to the "branched-polymer 
phase" (top), and the "crumpled phase" (bottom). 



Conversely, the bottom figure of Fig. 12 shows the effect of varying Z\, 
while keeping kq = 2.2 fixed. As A is decreased towards 0, we eventually hit 
another phase transition, separating the physical phase of extended universes 
from the CDT-equivalent of the crumpled phase of Euclidean quantum gravity, 
where the entire universe will be concentrated within a few time steps, as 
already mentioned above. (The point closest to the transition where we have 
taken measurements is the one in the bottom left-hand corner.) Also when 
approaching this phase transition the effective coupling constant fci goes to 
0, leading to the tentative conclusion that ki ^ along the entire phase 
boundary. 
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However, to extract the coupling constant G from (41) we not only have 
to take into account the change in ki , but also that in sq (the width of the 
distribution N^{i)) and in the effective four- volume C4, as a function of the 
bare coupling constants. Combining these changes, we arrive at a slightly 
different picture. Approaching the boundary where spacetime collapses in time 
direction (by lowering A), the gravitational coupling constant G decreases, 
despite the fact that l/fci increases. This is a consequence of §0 decreasing 
considerably. On the other hand, when (by increasing kq) we approach the 
region where the universe breaks up into several independent components, the 
effective gravitational coupling constant G increases, more or less like 
where the behaviour of ki is shown in Fig. 12 (top). This implies that the 
Planck length £pi = VG increases from approximately 0.48a to 0.83a when 
Ko changes from 2.2 to 3.6. Most likely we can make it even bigger in terms 
of Planck units by moving closer to the phase boundary. 

On the basis of these arguments, it seems likely that the non-perturbative 
CDT-formulation of quantum gravity does allow us to penetrate into the sub- 
Planckian regime and probe the physics there explicitly. Work in this direction 
is currently ongoing. One interesting issue under investigation is whether and 
to what extent the simple minisuperspace description remains valid as we go to 
shorter scales. We have already seen deviations from classicality at short scales 
when measuring the spectral dimension [25, 8], and one would expect them to 
be related to additional terms in the effective action (25) and/or a nontrivial 
scaling behaviour of fci . This raises the interesting possibility of being able to 
test explicitly the scaling violations of G predicted by renormalization group 
methods in the context of asymptotic safety [2]. 

4 Two-dimensional Euclidean quantum gravity 

The results described above are of course interesting and suggest that there 
might exist a field theory of quantum gravity in four dimensions (three space 
and one time dimension). However, the results are based on numerical sim- 
ulations. As already mentioned it is of great conceptional interest that we 
have a toy model, two-dimensional quantum gravity, where both the lattice 
theory and the continuum quantum gravity theory can be solved analytically 
and agree. Of course we can still be in the situation that there exists no de- 
scription of quantum gravity as a field theory in four dimensions (although 
we have presented some evidence in favour of such a scenario above), but we 
can then not blame the underlying formalism for being inadequate. 

4.1 Continuum formulation 

Let denote a closed, compact, connected and orientable surface of genus 

h and Euler characteristic x(^) = 2 — 2h. The partition function of two- 
dimensional Euclidean quantum gravity is formally given by 
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Z{A,G) = f2 [■D[g]e-''^''''^''\ (44) 

where A denotes the cosmological constant, G is the gravitational coupling 
constant and S is the continuum Einstein-Hilbert action defined by 



Sig; A,G)=A [ d^^y^ I d'^^R. 



(45) 



In eq. (44), we take the sum to include all possible topologies of two- 
dimensional manifolds (i.e. over all genera h), and in eq. (45) R denotes the 
scalar curvature of the metric g on the manifold M^. The functional integra- 
tion is over all diffeomorphism equivalence classes [g] of metrics on Al'^. 

In two dimensions the curvature part of the Einstein-Hilbert action is a 
topological invariant according to the Gauss-Bonnet theorem, which allows 
us to write 

oo 

Z(AG) = ^e'^W/«Z^(.l), (46) 

where 

Zh{A) = J V[g]e-^^^'^\ (47) 

and 

S{g;A) = AVg, (48) 

where Vg = f rf^^ ^ is the volume of the universe for a given diffeomorphism 
class of metrics. In the remainder of this section we will, for simplicity, restrict 
our attention to manifolds homeomorphic to S"^ or with a fixed number of 
holes unless explicitly stated otherwise. In this case we disregard the topolog- 
ical term in the action since it is a constant. The sphere S*^ with b boundary 
components will be denoted and we denote the partition function for the 
sphere, Zo{A) in eq. (47), by Z{A). 

In the presence of a boundary it is natural to add to the action a boundary 
term 

b 

S{g; A, Z„ Z^) = AVg +^ Z^L^^g, (49) 

i=l 

where Li^g denotes the length of the ith boundary component with respect 
to the metric g. We refer to the Zi's as the cosmological constants of the 
boundary components. The partition function is in this case given by 

WiA; Zu Z,) = J V[g] e-S(9;A,z„...,z,) _ ^^q) 

Since the lengths of the boundary components are invariant under diffeo- 
morphisms, it makes sense to fix them to values Li,...,Lf) and define the 
Hartle-Hawking wave junctionals by 
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W{A; Li, L,) = / V[g] e-^^^'^^ [] 5{Li - X,,,), (51) 

•' i=l 

where S{g;A) is given by eq. (48). Since eq. (50) is the Laplace transform of 
eq. (51), i.e. 

/•oo b 

W{A;Zu...,Zt)= / l[dLie-^''''W{A;L,,...,Lt), (52) 

•^0 i=l 

we denote them by the same symbol. We distinguish between the two by the 
names of the arguments. 



4.2 The lattice regularization 

At the outset wc restrict the topology of surfaces to be that of S'^ with a 
fixed number of holes. We view abstract triangulations of as defining a 
grid in the space of diffeomorphism equivalence classes of metrics on S^. Each 
triangle is a "building block" with side lengths a. This a will be a UV cut-off 
which we will relate to the bare coupling constants on the lattice. However, 
presently it is convenient to view a as being 1 (length unit). 

Let T denote a triangulation of 5^. The regularized theory of gravity will 
be defined by replacing the action Sg{A, Zi, . . . , Z\^ in eq. (49) by 

h 

5T(M,Ai,...,Afa) =MA^t + 5I^»'i' (53) 

1=1 

where Nt denotes the number of triangles in T and li is the number of links 
in the ith boundary component. The parameter /x is the bare cosmological 
constant and the Ai's are the bare cosmological constants of the boundary 
components. The integration over diffeomorphism equivalence classes of met- 
rics in eq. (50) becomes a summation over non-isomorphic triangulations. 
We define the loop functions (discretized versions of W{A, Zi, . . . , Zb j) by 
summing over all triangulations of S^: 

w{fi,X,,...,Xb)= E e-^-('''^^--'^^). (54) 

ii,-,ib Ter{h,...,h) 

Analogously, we define the partition function for closed surfaces by 

Z{u.)=Y^e-''-^'^\ (55) 

where Ct is the symmetry factor of T and SrifJ-) = fJ-Nt. Since we consider 
surfaces of a fixed topology wc have left out the curvature term in the action. 
It will be introduced later, when the restriction on topology is lifted. 
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Next we write down the regularized version of the Hartle-Hawking wave 
functional W{A, Li, . . . , Lf,): 

w{i^,h,...,h)= E ^"^"^"^ (56) 
TeT{h,...,it) 

with an abuse of notation similar to the one in the previous section. This can 
also be written in the form 

w{fi, h,...,lb) = Y, e~'"''wk,h,...,i, , (57) 

k 

where we have introduced the notation 

Wk,h,...,h 

for the number of triangulations in T{li. .... Ih) with k triangles. 

The discretized analogues of the Laplace transformations which relate 
W{A, Zi,...,Zb) and W{A, Li, . . . , Lb) are 

w{iJ,,Xi,...,Xb) = ^ e-^-^'^'^w{iJ,,h,...,lb) (58) 

and eq. (57). Similarly, we have for the partition functions 

Z{n)=J2e-'"'Z{k), (59) 

k 

Z{k)= E (60) 

TeT,Nt=k 

It follows from the definitions (57)-(58) that w{^, Ai, . . . , Ab) is the generating 
function for the numbers Wfe,;i,. the arguments of the generating function 
being and e~^'. In this way the evaluation of the loop functions of two- 
dimensional quantum gravity is reduced to the purely combinatorial problem of 
finding the numher of non-isomorphic triangulations of S'^ or S'^ with a given 
number of triangles and boundary components of given lengths. 
We use the notation 

w{g,zu...,Zb)= Yl «^Mi,...^/^r''"^---^6"'''"^ (61) 

for the generating function with an extra factor z^^.-.z^^, i.e. we make the 
identifications 

g = e-^, Zi = e^\ (62) 

The reason for this particular choice of variables in the generating function is 

motivated by its analytic structure, which will be revealed below. 

In the following we consider a particular class of triangulations which in- 
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Fig. 13. A typical unrestricted "triangulation" . 



eludes degenerate boundaries. It may be defined as the class of complexes 
homeomorphic to the sphere with a number of holes that one obtains by suc- 
cessively gluing together a collection of triangles and a collection of double- 
links which we consider as (infinitesimally narrow) strips, where links, as well 
as triangles, can be glued onto the boundary of a complex both at vertices 
and along links. Gluing a double link along a link makes no change in the 
complex. An example of such a complex is shown in Fig. 13. The reason we 
use this class of triangulations is that they match the "triangulations" we 
obtain from the so-called matrix models to be considered below. We call this 
class of complexes "unrestricted triangulations" . 

One could have chosen a more regular class of triangulations, correspond- 
ing more closely to our intuitive notion of a surface. However the degenerate 
structures present in the unrestricted triangulations appear on a slightly larger 
scale in the regular triangulations in the form of narrow strips consisting of 
triangles. Since wc want to take the lattice side a of a triangle to zero in the 
continuum limit, there should be no difference in that limit between various 
classes of triangulations, unless more severe constraints are introduced. We 
say that "the continuum limit is universal" . But at some point the constraint 
can be so strong that the continuum limit is changed. We will meet precisely 
such a change below, leading from (Euclidean) dynamical triangulations (DT) 
to causal dynamical triangulations (CDT). 

Let w{g, z) denote the generating function for the (unrestricted) triangu- 
lations with one boundary component. Then we have 



We have included the triangulation consisting of one point. It gives rise to 

the term 1/z and wc have tii{){g) = 1. The function wi{g) starts with the 
term g, which corresponds to an unrestricted triangulation with a boundary 




(63) 



k=0 1=0 



1=0 
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consisting of one (closed) link with one vertex and containing one triangle. 
The coefficients w^ i in the expansion Wi^ig + ig^ + • • • of ^1(5) are the 
numbers of unrestricted triangulations with a boundary consisting of one link. 

The coefficients w{g,z) fulfill a recursion relation which has the simple 
graphical representation shown in Fig. 14. The diagrams indicate two opera- 




Fig. 14. Graphical representation of relation (64): The boundary contains one 
marked link which is part of a triangle or a double link. Associated to each tri- 
angle is a weight g, and to each double link a weight 1. 



tions that one can perform on a marked link on the boundary to produce a 
triangulation which has either fewer triangles or fewer boundary links. The 
first term on the right-hand side of Fig. 14 corresponds to the removal of a tri- 
angle. The second term corresponds to the removal of a double-link. Note that 
removing a triangle creates a new double-link if the triangle has two boundary 
links. In addition, note that we count triangulations with one marked link on 
each boundary com,ponent and adopt the notation introduced above for the 
corresponding quantities. 

The equation associated with the diagrams is 



M9,z)]k, = [gzw{g,z)]^:, + 



-w'^{g,z) 



(64) 

k,l 



The subscripts k, I indicate the coefficient of g^ jz^^^ . Let us explain the equa- 
tion in some detail. The factor gz in eq. (64) is present since the triangulation 
corresponding to the first term on the right-hand side of Fig. 14 has one tri- 
angle less and one boundary link more than the triangidation on the left-hand 
side. The function w'^{g.z) in the last term in eq. (64) arises from the two 
blobs connected by the double-link in Fig. 14 and the 1/z in front of w'^{g, z) 
is inserted to make up for the decrease by two in the length of the boundary 
when removing the double link. 

As the reader may have discovered, eqs. (64) is not correct for the smallest 
values of I. Consider Fig. 14. The first term on the left-hand side of eq. (64) (a 
single vertex) has no representation on the diagram. In order for eq. (64) to be 
valid for A; = ^ = we have to add the term l/z on the right-hand side of eq. 
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(64). Furthermore, it is clear from Fig. 14 that the first term on the right-hand 
side has at least two boundary links. Consequently, the term gzw{g, z) on the 
right-hand side of eq. (64) should be replaced by gz{w{g, z) — 1/ z — Wi{g)/ z"^) 
such that all terms corresponding to triangulations with boundaries of length 
and 1 are subtracted. It follows that the correct equation is 



We will refer to eq. (65) as the loop equation. It is a second-order equation 
in w((7, z). As will be clear in the following this algebraic feature allows us to 
extract asymptotic formulas for the number of triangulations with k triangles 
in the limit fc — > oo. 

4.3 Counting graphs 

Let us begin by solving eq. (65) in the limit g = Q. In this case there are 
no internal triangles and the triangulations arc in one-to-one correspondence 
with rooted branched polymers ^ . The double-links correspond to the links of 
the branched polymers and the root is the marked link, see Fig. 15. If = 



Fig. 15. Rooted branched polymers created by gluing of a boundary with one 
marked link. 



® One might think that such polymers arc not relevant at all for studying real sur- 
faces made of triangles, not to mention higher pieccwisc linear manifolds, but in 
fact the branched polymer structure is quite generic. Surfaces or higher dimen- 
sional manifolds can "pinch", such that two parts of the triangulation are only 
connected by a minimal "neck" . If this happens in many places one can effectively 
obtain a branched polymer structure even for higher dimensional picccwise linear 
manifolds. Such minimal necks have been used to measure critical exponents of 
various ensembles of piecewise linear manifolds [27] and in four-dimensional Eu- 
clidean quantum gravity one has indeed, as mentioned above, observed a phase 
where the four-dimensional piecewise linear manifolds degenerate to branched 
polymers [23] . The same is the case for bosonic strings with central charge c > 1 



{z- gz'^)w{g,z 



')-! + g{wi{g) + z) = w^{g, z). 



(65) 




[28]. 
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then eq. (65) reads 

w'^{z) - zwiz) + 1 = 0. (66) 

The above equation has two sohitions. The one that corresponds to the eount- 
ing problem has a Taylor expansion in whose first term is z^^ (recall that 
Wofi = 1). This solution is given by 

wiz) = l{z- V^^) . (67) 
Expanding in powers ofl/z yields 

oo 

1=0 ^ 

where 

and W21 is the number of rooted polymers with I links. Note that W21 are the 
Catalan numbers, known from many combinatorial problems. 

The generating function w{z) is analytic in the complex plane C with a cut 
on the real axis along the interval [—2, 2]. The endpoints of the cut determine 
the radius of convergence of w{z) as a function of 1/z or, equivalently, the 
exponential growth of W21 ■ 

We can solve the second-order equation (65) and obtain 

wig, z) = l (v'iz) - ViV'izW - ^Q{z)) , (70) 

where, anticipating generalizations, we have introduced the notation 

V'{z) = z-gz'', Q{z) = l-gwi{g)-gz. (71) 

The sign of the square root is determined as in eq. (67) by the requirement 
that w{g,z) = 1/z + 0(1/ z^) for large z (since wq.o = 1). If g = then 
V'{zf - 4Q{z) = z'^ -A. For 5 > 0, on the other hand, V'{zf - 4Q{z) is a 
fourth-order polynomial of the form 

V'izf - 4Q{z) = {z-{2 + 2g) + 0{g^)} 

x{z+{2- 2g) + 0{g^)}{gz - (1 - 2g^) + 0{g^)Y (72) 

in a neighbourhood oi g = Q since the analytic structure of w{g,z) as a 
function of z cannot change discontinuously at = 0. We can therefore write 

V'izf - AQ{z) = {z- c+{g)){z - c.{g)){c2{g) - gzf, (73) 

and, by eq. (70), 
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w{9, = 9z^ + [gz - C2)\/ {z - c^){z - c_)) , (74) 

where c_, c+ and C2 arc functions of g, analytic in a neighbourhood of g = 0. 
We label the roots so that c_ < c+. The numbers c_, c+ one? C2 are uniquely 
determined by the requirement that w{g, z) = 1/ z + 0{l/ z"^), again originating 
from, ?y'o.o = 1- This requirement gives three equations for the coefficients of 

We can generalize the above counting problem to planar complexes made 
up of polygons with an arbitrary number j < n of sides, including "one-sided" 
and "two-sided" polygons. If we attribute a weight gtj to each j-sided polygon 
and a weight z to each boundary link, and adopt the notation 

V'{z) = z- g{ti + t2Z + taz'^ + ■■■+ tnz""-^), (75) 

n 0-2 

Q{z) = 1 - ^t,- ^2'«;,_2-K5), Ma) = 1, (76) 

j=2 1=0 

the analogue of eq. (65) is 

w{g,zf = V'{z)w{g,z)-Q{z), (77) 

or 

w{g, z) = g ^i^+t2+t3Z-\ h in-j""^^ w{g, z) 

+-Q{z) + -w\g,z). (78) 
z z 

The graphical representation of (78) is shown in Fig. 16. The subtraction of 
the polynomial Q{z) in cq. (77) reflects the fact that the term with a j-sidcd 
polygon in Fig. 16 must have a boundary of length at least j — 1 for j > 1. 
The constant term 1 in Q corresponds to the complex consisting of a single 
vertex. 

The solution can be written 

w{9,z) = \ (y'{z) - M{zW{z-c+{g)){z-c.{g)) , (79) 

where M{z) is a polynomial of a degree which is one less than that of V'{z). 
Again, the polynomial M is uniquely determined by the requirement that 
w{g,z) falls off at infinity as before, i.e. w{g,z) = 1/z -\- 0{l/z^), and the 
additional requirement that 'w{z) has a single cut. It is sometimes convenient 
to write (79) as 

w{g, z) = \ {v'{z) - X5 Mk{g){z - c+f-' ^{z-c+{g)){z-c-{g)^ . 

(80) 
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Fig. 16. Graphical representation of relation (78): the marked link of the boundary 
belongs either to an i-gon (associated weight g ti) or a double link (associated weight 
1). It is also a graphical representation of eq. (114) if instead of weight 1 we associate 
a weight ga to the marked double link. 



One can show the following: for ti, . . . , > and f„ > one has 

Mk{g) < 0, k<n, (81) 

while 

Ml (5) > (82) 

in a neighbourhood of 5 = 0. When we increase g we first reach a point Qc, 
where (81) is still satisfied but 

Mi(ffe) = 0. (83) 

The coupling constant point Qc is thus the point where the analytical structure 

of w{g, z) changes from being identical to that of the branched polymer, i.e. 
it behaves like {z - c+Y^'^iz - c-)l/^ to {z - c+f/'^{z - c-fl"^. The function 
w{g^ is an analytic function around the point 9 = {w{g = 0, z) is the 
branched-polymer partition function discussed above). The radius of conver- 
gence is precisely gc- If we return to the expansion in (64) each term wi{(j) 
has this radius of convergence. It is the generating function for triangulations 
with one boundary consisting of I links. The singularity of wi{g) for g — > g^, 
determines the asymptotic behaviour of the number of such triangulations for 
a large number of triangles, i.e. the leading behaviour of the numbers w^^^i for 
large k (see (89) and (90) below). 

Let us introduce g^ = gtj as new variables. We have 

n 

w{guz)= Yl Hk,},iz-^'+'^l[g-\ (84) 

l,ki,...,kn j = l 

where W{kj},i is the number of planar graphs with kj j-sided polygons, j = 
1, . . . , n, and a boundary of length 
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Prom w{gi,z) we can derive the generating function for planar graphs with 
two boundary components by applying the loop insertion operator 

dV{z)=^^J^dij- ^^^^ 

One should think of this operator as acting on formal power series in an ar- 
bitrary number of variables gj. The action of d/dV{z2) on w{gi,Zi) has in 
each term of the power series the effect of reducing the power kj of a specific 
coupling constant gj by one and adding a factor jkj/z^^^. The geometrical 
interpretation is that a j-sided polygon is removed, leaving a marked bound- 
ary of length j to which the new indeterminate Z2 is associated. The factor 
kj is due to the possibility to make the replacement at any of the kj j-sided 
polygons present in the planar graph, while j is the number of possibilities 
to choose the marked link on the new boundary component. The generat- 
ing function for planar graphs with b boundary components can therefore be 
expressed as 

w{9,, .1, ...,.,) = ^ ■ ■ ■ ^-(5., z,). (86) 

A most remarkable result is the following: for any potential V{gi,z) the 
two-loop function w{gi, Z\,Z2) has the form 

w{gi,zi,Z2) (87) 

_ 1 / ZiZ2 - ^{Zl+ Z2){c+ + C-) + C+C- \ 

- 2{Z, - Z2f \^[{Z, - C+)(Z1 - C.)][{Z2 - C+){Z2 - )' 

Note that there is no explicit reference to the potential V{gi, z), but of course 
c+ and c_ depend on the potential. 

Prom this formula one can in principle construct the multi-loop function 
'w{g, zi, . . . ,Zb) by applying the loop insertion operator 6—2 times. One can 
use this formula to find the leading singularity w{g, Zi, . . . , Zb) when g gc, 
the critical value of the coupling constant g and the value where Mi{g) = 0. 
One finds 

^ \ 26-5 



K,,.„...,.,)~^^_j , (88) 

as 5 — > gc- This implies that the generating function w{g,li, . . . ,lb) for 
the number of triangulations, Wk^ii,...,i,,, constructed from k triangles with 
b boundary components of length li, . . . ,lb, has a singularity as g ^ gc that 
is independent of the length of the boundary components and is given by 

^ X 26-5 



Whig,li,...,lb)^ (^^==^] . (89) 
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Finally, we obtain from eq. (89) the asymptotic behaviour of Wfe,;!,...,;^ as 
fc ^ oo: ^ 

WkM,-,h ~ (^) (90) 

We note that these results can be generalized to triangulations which have h 
handles: 

^ ^ 26+(/(-l)5 



Wh{g,li,...,h)-^ ( ^/g-r^ l • (91) 



- ( I ) fc^^'-^^^^^-^ (92) 



For future applications it is important to note that the position of the leading 
singularity gc in (91) or, alternatively, the exponential growth of the number 
of triangles in (92), is independent of the number of handles or the number of 
boundaries. 



4.4 The continuum limit 

We now show how continuum physics is related to the asymptotic behaviour 
of w'k}i,...,h for A: ^ 00 and ii, . . . , 4 ^ c)o in a specific way, and we use the 
results for the generating functions Wh{g,z\, . . . , Zb) derived in the previous 
sections to study this limit. 

Before discussing details it is useful to clarify how we expect the contin- 
uum wave functional W{A, Zi, . . . , Zb) to renormalize. Since the cosmological 
constants A and Zi have dimensions and 1/a, respectively, a being the 
length of the lattice cutoff, it is natural to expect that they are subject to an 
additive renormalization 

Ae = ^+A, Zi^, = hhl+Zi, (93) 

where Ac and Zi^c are the bare cosmological coupling constants. Since our 
regularization is represented in terms of discretized two-dimensional mani- 
folds, the bare cosmological constants should be related to the dimensionless 
coupling constants fi, A, by 

^=4, Z,,, = ^, (94) 

so that eq. (93) can be written 

jjL- lic = a^A, Xi - Xi^c = aZi. (95) 

In the following we assume for simplicity that all the Aj,cS are equal to Ac. We 
identify the constants /Uc and Ac with the critical couplings Qc and c+(fifc) via 
the relations 
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-=e'^% c+{g,) = e^'. (96) 
9c 

Recalling the relation (62) between jjL, g and z, A, it follows that the a ^ limit 
of the functions w{pL,\i, . . . ,\h) is determined by their singular behaviour 
at gc- The renormalization (93) has the effect of cancelling the exponential 
entropy factor for the triangulations, see eq. (90). Note that since we have the 
same exponential factors for all genera, we expect the renormalization of the 
cosmological constants to be independent of genus. 

We begin by studying the continuum limit for planar surfaces. Then we 
will discuss how to take higher genera into account, thereby reintroducing 
the gravitational coupling constant G and also discussing its renormalization. 
This will lead us to the so-called double-scaling limit. 

We are interested in a limit of the discretized models where the length a 
of the links goes to zero while the number k of triangles and the lengths Zj of 
the boundary components go to infinity in such a way that 

V = ka^ and = ka (97) 

remain finite. The asymptotic behaviour of Wk^i-^ is given by cq. (90) if 

the h,. . . ,lb remain bounded. In this case the leading term is of the form 

where /? is a critical exponent. If the boundary lengths diverge ac- 

cording to eq. (97), we expect a corresponding factor 

e^«-^^/"(Li/a)", 

where a is another critical exponent. This form of the entropy was encountered 
for branched polymers in eq. (69). We can therefore express the expected 
asymptotic behaviour of the coefficients 'Wk,ii,...,ii, as 



.^v E. h a-"b-2/3 ^r{y, Li,...,Lh), (98) 



as a ^ 0, with V and Li defined by eq. (97) fixed. The factor a~°'^~'^^ may 
be thought of as a wave-function renormalization. 

From cq. (98) we deduce that the scaling behaviour of the discretized wave 
functional ■u;(/U, Ai, . . . ., A^) is given by 

«;(/x,Ai,...,A6) = e-'^'=e-^^^'''Wfe,i,, 

k,li ,...,lb 



_1_ ^ e-('^-'^°)'=e-^.(^'-^°)'-VF(F,Li,...,Lb) 

W{A,Zi,...,Zb), (99) 



1 

g(a+l)fc+(2/3+2) 

where we have used the relation 
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W{A,Z^,...,Z,)= [ dl/TTdLie-^^-^.^'^- W^(y,Li,...,X6). (100) 

Our next goal is to show that we can take a hmit as suggested by eqs. (95) 
and (99). In terms of the variables g, Zi we have 



9 = 50(1-^0^), Zi = Zc{l + aZi), 
where we have introduced the notation 

Zc = c+{gc) = e^" 



(101) 



(102) 



for the critical value c+(gc) of z corresponding to the largest allowed value of 
g. Inserting (101) and (102) in the expression (87) one obtains 



where 



W{A,Z^,Z2) 



w{g,zi,Z2)r^a-^W{A,Zi,Z2) 



+ Z2) + VA 



^ - ^2)' \yJ{Z, + VA){Z2 + VA) 



(103) 



(104) 



Similarly one can show, using the loop insertion operator, that when the 
number of boundaries is larger than two one has 



w{g,zi,...,Zb) 



1 



76/2-5 



a 



I.e. 



W{A,Z,,...,Zb) 



d_ 
"dA 



d_ 
dA 



&-3 



b-3 



y/AfJi{Zi + VAf/^ 



(105) 



1 

7i 



n 



\ {Zi + x/I)3/2 



The continuum expressions for the n-loop functions are all universal and in- 
dependent of the explicit form of the potential V{z) as long as the weights 
U > 0. For the one-loop function the situation is different. As is seen by 
formally applying the counting of powers a in (105) to the case 6=1 one 
obtains the power a^/^, i.e. a positive power of a. The important point in 
(105) and (103) is that the power is negative: in the scaling limit a ^ these 
terms will dominate. This is how the formulas should be understood: there 
arc other terms too, but they will be subdominant when a 0, i.e. when 
g g,; and z ^ Zc &s dictated by (101) and (102). For the one- loop function 
the term associated with a'^/^ will vanish when a ^ and wo will be left with 
a non-universal term explicitly dependent on the potential V. However, the 
term associated with a^/^ is still the leading term which is non-analytic in the 
coupling constant g, so if we differentiate a number of times with respec;t to 
g and then take the limit a ^ it will be dominant. No continuum physics 
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is associated with the analytic terms since they contain g only to some finite 
positive power, and arc thus associated with only a finite number of triangles 
(of which the lattice length a ^ when g g^) if we recall the interpretation 
of w{g, z) as the generating function of the number of triangulations. Inserting 
(101) and (102) in the expression (80) one obtains 

w{g, z) = ^iV'iz) + a^/^W{A, Z) + 0(aV2)) (106) 
where 



W{A, Z) = {Z- ^^)\fzWI. (107) 

This ends the calculation of the loop-loop correlation functions for mani- 
folds with topology 5^, the sphere with b boundaries. The results agree with 
continuum calculations using quantum Liouville theory and it follows that one 
can obtain a continuum, diffeomorphism-invariant theory starting out with a 
suitable lattice theory, where the lattice link length acts as a diffeomorphism- 
invariant UV cut-off and simply taking the lattice spacing a — > 0, while renor- 
malizing the bare couplings in a standard way, namely, the cosmological and 
the boundary cosmological coupling constants. 

Let us end this description of Euclidean quantum gravity by mentioning 
the corresponding results for higher-genus surfaces. The generalization of (105) 
is ^ 

Wh{A,Z,,...,Zb). (108) 
In particular, taking 6 = leads to the expression 

Z'^ia) - (^.^ip^, (109) 

for the singular part of the partition function, where the constants th can in 
principle be (and have been) computed. 

In (109) we have actually completed the task of calculating (47). We can 
now reintroduce the gravitational coupling constant G and try to calculate 
the complete partition function 

•oo 

X ^ 2 — 2h r',^ I \ 5(1 — h) 

Z{G,A) = J2The—a^^^-''U^^. (110) 

The factor present for each genus can be absorbed in a renormalization 
of the gravitational coupling constant 

1 1 5, yl , , 

-7log-o, (111) 



Gren G{a) 

whore Gren dcnotc;s the renormalized gravitational coupling. A continuum 
limit of (110) only exists in the limit a ^ if we allow G to be a function of 
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the lattice spacing a determined by (111) for fixed Gj-en and A. The continuum 
limit is then given by 

Z{G,A)=Y,rh [^-^""A-i) , (112) 

and depends only on the variable x — Ae^/'-^'^\ To actually calculate Z{G, A) 
we have to perform the summation over the number of handles h in (112), an 
interesting task which we will not address here. Rather, we will focus on eq. 
(Ill), since we can use this equation to calculate the /3-function for G using 

= - 52G^. (113) 

Two-dimensional Euclidean quantum gravity is asymptotically free as already 
mentioned in the introduction. 



da 



5 Two-dimensional Lorentzian quantum gravity 

As already mentioned above, Euclidean quantum gravity does not really work 
in more than two dimensions (in the sense of leading to a continuum theory of 
higher-dimensional geometry). By contrast, the formalism called CDT, based 
on causal dynamical triangulations, seems to lead to very interesting results. It 
is based on the idea that there exists a globally defined (propor-)time variable, 
which can bo used to describe the evolution of the universe. In addition, one 
assumes that the topology of space is unchanged with respect to the foliation 
defined by this global time. 

These requirements are definitely not satisfied in two-dimensional Eu- 
clidean quantum gravity. In principle one can also "superimpose" a proper 
time on Euclidean quantum universes and follow their evolution as first de- 
scribed in the seminal work by Kawai and collaborators. Starting out with a 
spatial universe of topology S^, it will immediately split up into many dis- 
connected spatial , one-dimensional universes as a function of proper time. It 
turns out that the structure is fractal, in the sense that an infinity of spatial 
universes, most of them of infinitesimal spatial extension, will be created as a 
function of proper time. 

Since two-dimensional Euclidean quantum gravity is explicitly solvable, 
even on a lattice before the continuum limit is taken, as described above, it 
is of interest to understand the transition from the Euclidean lattice gravity 
theory to the CDT lattice gravity theory. Clearly one has to suppress the 
splitting of a spatial universe into two or more disconnected spatial universes 
if one wants to move from the spacetime configurations which characterize 
the Euclidean path integral to the configurations present in the CDT path 
integral. It makes sense to talk about the splitting of a spatial universe into 



46 



J. Ambj0rn, J. Jurkiewicz and R. Loll 



two if the universe has Lorentzian signature, since such a spUtting (in the 
simplest case) is associated with an isolated point where the metric and its 
associated light-cone structure are degenerate, which has a diffeomorphism- 
invariant meaning. This was the motivation for imposing such a constraint 
in the original CDT model. By working in Lorentzian signature initially, and 
only later rotating to Euclidean signature, this constraint survives also in 
(the Euclidean version of) CDT. Going back to Fig. 16, this suggests that 
one should associate a factor gs instead of a factor 1 with the graph with 
the double line. Geometrically this figure can be viewed as a process where 
a triangle is removed at a marked link (and a new link is marked at the 
new boundary), except in the case where the marked link does not belong 
to a triangle, but is part of a double-link, in which case the double link is 
removed and the triangulation is separated into two. If one thinks of the 
recursion process in Fig. 16 as a "peeling away" of the triangulation as proper 
time advances, the presence of a double link represents the "acausal" splitting 
point beyond which the triangulation splits into two discs with two separate 
boundary components (i.e. two separate one-dimensional spatial universes). 
The interpretation of this process, advocated in [33], is that it represents a split 
of the spatial boundary with respect to (Euclidean) proper time. Associating 
an explicit weight with this situation and letting gs suppress this 
process compared to processes where we simply remove an i-gon from the 
triangulation. Nevertheless, we will see below that there exists an interesting 
scaling of gs with a such that the process survives when we let a ^ 0, but 
with a result different from the Euclidean quantum gravity theory. We call 
this new limit generalized CDT [17]. 

Let us introduce the new coupling constant Qs in eq. (78). The equation is 
then changed to 



In the analysis it will be convenient to keep the coupling constant ti > 0, 
although we are usually not so interested in situations with one-gons. It can 
be motivated as follows. Consider a "triangulation" consisting of Ti one-gons, 
T2 two-gons, T3 triangles, T4 squares etc., up to T„ n-gons. The total coupling- 
constant factor associated with the triangulation is given by 



We observe that in the limit gs ^ 0, a, necessary condition for obtaining a 
finite critical value gdds) is Ti > 0. We should emphasize that the analysis 
described below can be carried out also if we suppress the appearance of any 
one-gons (by setting ti = 0) in our triangulations, but it is slightly more 
cumbersome since then gdds) — > 00 as — * 0, requiring further rescalings. 




(114) 



,Ti + ---+T„ -Ti/2-|-T3/2+--- + (n/2-l)T„ 



(115) 
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For simplicity we will consider the simplest nontrivial model with poten- 
tial*^ 

and analyze its behaviour in the limit Qs — > 0. The disk amplitude (79) now 
has the form 

w{z) = ^(^-g + z-gz^+ g{z - C2W{z - c+){z - c_)) , (117) 

and the constants C2, c+ and c_ are determined by the requirement that 
w{z) — > 1/2; for 2; — > oo. Compared with the analysis of the previous section, 
the algebraic condition fixing the coefficient of 1/z to be unity will now enforce 
a completely different scaling behaviour as ^ 0. 

For the time being, we will think of Qs as small and fixed, and perform 
the scaling analysis for gdgs)- As already mentioned above, the critical point 
gc is determined by the additional requirement Mi = in the representation 
(80), i.e. that 02(^0) = c+{9c), which presently leads to the equation 

{l-4g^f^' = uV3glgs. (118) 

Anticipating that we will be interested in the limit gs 0, we write the 
critical points as 

9c{9s) = i(i - ^gdgs)), ^9c{gs) = IgV^ + 0(9^/% (ii9) 

and 

(120) 

while the size of the cut in (116), c+(gc) — c-{gc), behaves as 

c+(5e) - c_(5c) = 45]/^ + 0{g'J'). (121) 

Thus the cut shrinks to zero as 9,, 0. 

Expanding around the critical point given by (119)-(120) a nontrivial limit 
can be obtained if we insist that in the limit a — > 0, gfg scales according to 

gs = G,a^ (122) 

where a is the lattice cut-off introduced earlier. With this scaling the size of 
the cut scales to zero as AaCV^. In addition \J{z — c+)(2 — c_) a a if we 
introduce the standard identification (101): z = c+(gc) + aZ. This scaling is 
different from the conventional scaling in Euclidean quantum gravity where 

® the rationale for calling V a "potential" will become clear below 
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\J{z — C-|-)(2; — c_) oc a}-!'^ since in that case (z — c+) scales while (z — cJ) 
does not scale. 

We can now write 

9 = 5c(ff«)(l - a^A) = g{l - a^^^dt + 0(a')), (123) 
with the identifications 

A,^,^A+^Gy\ g=^, (124) 

as well as 

z = Zc + aZ = z + aZ^di + Oia^), (125) 
with the identifications 

Z,dt = Z + Gy\ z = l. (126) 
Using these definitions one computes in the limit a ^ that 



In (127), the constant H (or rather, its rescaled version h = H/\/2A^ ) 
satisfies the third-order equation 

which follows from the consistency equations for the constants C2, c+ and c_ 
in the limit a — > 0. We thus define 

wiz) - - Wcdt(^cdt,^dt,G,) = -WiZ,A,Gs) (129) 

a. a 

in terms of the continuum Hartle-Hawking wave functions Wcdt(-^cdt5 ^cdt, Gg) 

and W(Z, A, G^). 

Notice that while the cut of y/{z ~ c+){z — C-) goes to zero as the lattice 
spacing n. it nevertheless survives in the scaling limit when expressed in terms 
of renormalized "continuum" variables, as is clear from eq. (127). Only in the 
limit Gg ^ it disappears and we have 

«;(2) = Jw^cdt(^cdt,^cdt,G,) — ^ 7 J/OA— ' (130) 

a Ga^o a Zcdt + V ^^cdt 

which is the original CDT disk amplitude introduced in [10]. 
Let us make some comments: 

(1) We have dealt here directly with a "generalized" CDT model which in 
the limit Gg ^ reproduces the "original" CDT disk amplitude (130). The 




original two-dimensional CDT model was defined according to the principles 
already outlined in our discussion of four-dimensional quantum gravity. Thus 
the interpolation between two spatial slices separated by one lattice spacing 
of proper time is shown in Fig. 17. This figure is the two-dimensional analogue 
of Fig. 4. The combinatorial problem of summing over all such surfaces con- 
necting two spatial boundaries separated by a certain number of time steps 
can be solved [10], and the continuum limit can be taken. The corresponding 
continuum "propagator" is found to be 

4/1 dte~^^^^ 
1 

(V^^ + x)(^/z;^+y)-e-2v:^r(yz;^-x)(x/z;^-y)' 

where X, Y arc the boundary cosmological constants associated with the two 
boundaries of the cylinder and T is the proper time separating the two bound- 
aries. The propagator has an asymmetry between X and Y because we have 
marked a point (a vertex in the discrctizcd model). By an inverse Laplace 
transform one can calculate the propagator G{X, L; T) as a function of the 
length of the unmarked boundary. In particular, we have 

/ioo 
dYG{X,Y-T), (132) 
-ioo 

and we define the CDT disk amplitude as 

pOO 

H^(o)(X)= / dT G{X,L = Q-T), (133) 
Jo 

and it is given by the expression on the far right in (130). 

The generalized CDT model allows branching of the spatial universes as 
a function of proper time T, the branching being controlled by the coupling 
constant Gs- This results in the graphical equation for the generalized CDT 
disk amplitude shown in Fig. 18. The corresponding equation is: 
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Fig. 18. Graphical illustration of eq. (134). Shaded parts represent the generalized 
CDT disk amplitude, unshaded parts the original CDT disk amplitude and the 
original CDT propagator (131). 



W{X) = W^^HX) + Gs dT dLidLa 
Jo Jo 

(Li + L2)GiX, Li + L2; T)W{L^)W{L2), 



(134) 



where W{L) is the disk amplitude corresponding to boundary length L. It can 
be solved [17] for W[X) and the solution is Wcdt(-'S^, ^cdt, Gs) found above. 

(2) Using eq. (128) we can expand w{z) into a power series in Gs/(2ylcdt)'^^^ 
whose radius of convergence is l/3-\/3. For fixed values of /Icdt, this value cor- 
responds to the largest value of Gs where (128) has a positive solution for h. 
The existence of such a bound on Gg for fixed ylcdt was already observed in 
[17]. This bound can be re-expressed more transparently in the present Eu- 
clidean context, where it is more natural to keep the "Euclidean" cosmological 
constant A fixed, rather than Acdt- We have 



G., 



< 



1 



3G: 



■2/3. 



(2yledt)3/2 - 3^ + 3G: 



2/3 



< 1, 



(135) 



which for fixed yl > is obviously satisfied for all positive Gs . In order to see 
that the usual Euclidean 2d quantum gravity (characterized by some finite 
value for gs) can be re-derived from the disk amplitude (127), let us expand 
eq. (127) for large Gg- The square root part becomes 



(136) 



which coincides with the generic expression a^^^W{Z, A) in Euclidean 2d 
quantum gravity (c.f. eqs. (106) and (107)) if we take Gs to infinity as gs/a^ 
(and take into account a trivial rescaling of the cosmological constant). How- 
ever, if we reintroduce the same scaling in the l/'(z)-part of 'w{z), it does not 
scale with a but simply goes to a constant. This term would dominate 'w{z) 
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in the limit a ^ if one did not remove it by hand, as is usually done in the 
Euclidean model. 

(3) Why does the potential V'{z) (and therefore the entire disk amplitude 
w{z)) scale (like 1/a) in the new continuum limit with = GgO?, a ^ 0, 
contrary to the situation in ordinary Euclidean quantum gravity? This is 
most clearly seen by looking again at the definitions (123) and (125). Because 
of the vanishing 

V'{z,g) = Q, V"{z,g)=Q (137) 

in the point [z^g) = (1, 1/2), expanding around {z,g) according to (123), (125) 
leads automatically to a potential which is of order a? when expressed in terms 
of the renormalized constants (.^cdt) ^cdt)) precisely like the square-root term 
when expressed in terms of (Zcdt,^cdt)- 

The point (z,.g) differs from the critical point (zc(.9s), .9c(.9s)), as long as 
£(s ^ 0. In fact, both 1 /z and g lie beyond the radii of convergence oil/ z and g, 
which are precisely l/zc{gs) and gdgs)- However, since the differences are of 
order a and , respectively, they simply amount to shifts in the renormalized 
variables, as made explicit in eqs. (124) and (126). Therefore, re-expressing 
W{Z, A, Gs) in (129) in terms of the variables Z^dt and ylcdt simply leads to the 
expression Wcdt('^cdt, ^cdt, G^), first derived in [17]. Similarly, any geometric 
quantities defined with respect to Z and A can equally well be expressed in 
terms of Zcdt and Acdt- For instance, the average continuum length of the 
boundary and the average continuum area of a triangulation are given by 

dlnW{Z,A,Gs) dlnW,dt{Zcdt,A,dt,Gs) 

= — dz — = dz:;, ' ^'^^^ 

_ dlnW{Z,A,Gs) _ dlnW,dt{Z,dt,A,dt,Gs) 
^ ^ ~ dA - 9^^^ • ^'"^^^ 

In the limit of Gs 0, the variables {Z,A) and (.^cdti^cdt) become iden- 
tical and the disk amplitude becomes the original CDT amplitude {Z^dt + 
a/21^)"^ alluded to in (130). 



6 Matrix model representation 

Above we have solved the two-dimensional gravity models by purely combina- 
torial techniques which emphasize the geometric interpretation: the quantum 
theory as a sum over geometries. The use of so-called matrix models allows one 
to perform the summation over the piecewise linear geometries in a relatively 
simple way. Surprisingly, it turns out that the scaling limit of the generalized 
CDT model has itself a matrix model representation. We will here describe 
how matrix models can be used instead of the combinatorial methods and 
how one is led to the CDT matrix model. 
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Fig. 19. The matrix representation of triangles which converts the gluing along 
links to a Wick contraction. 

Let be a Hermitian NxN matrix with matrix elements (pai3 and consider 
for fc = 0, 1, 2, . . . the integral 

y^e-i^'-^^l^itr^^y^ (140) 

where 

# = dRC(pal3 Y\_ dlm(l)al3. (141) 
a</3 a</3 

Wc can regard <p as a zcro-dimcnsional matrix-valued field so the integral can 
be evaluated in the standard way by doing all possible Wick contractions of 
(tr (f)^)'' and using 

{M«'0') = cJd<P e-^^"^ = Sa0'S0a', (142) 

where C is a normalization factor. The evaluation of the expression (140) 

can be interpreted graphically by associating to each factor tr cj)^ an oriented 
triangle and to each term 4>a04'i3-f4'-fa contributing to the trace a labeling of 
its vertices by a,/3,^ in cyclic order, such that the matrix element ^ai3 is 
associated with the oriented link whose endpoints are labeled by a and (3 in 
accordance with the orientation, eq. (142) can then be interpreted as a gluing 
of the link labeled by {a/3) to an oppositely oriented copy of the same link, 
see Fig. 19. 

In this way the integral (140) can be represented as a sum over closed, 
possibly disconnected, triangulations with k triangles. Triangulations with an 
arbitrary genus arise in this representation. By summing over k in (140) it fol- 
lows by standard arguments that the (formal) logarithm of the corresponding 
integral is represented as a sum over all closed and connected triangulations. 
The contribution of a given triangulation can be determined by observing that 
in the process of gluing we pick up a factor of N whenever a vertex becomes 
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an internal vertex in the triangulation. Thus, the weight of a triangulation T 
is simply 

If we make the substitution 

tr<^3__^tr<^3^ (143) 
the weight of T is replaced by 

gkjqN,{T)-k/2 ^ gkjqx(T) ^ 

where x{T) is the Euler characteristic of T. Note also that the factor (fc!)"""^ 
ill (140) is cancelled in the sum over different triangulations because of the 
k\ possible permutations of the triangles, except for triangulations with non- 
trivial automorphisms, in which case the symmetry factor C^^ survives. With 
the identifications ^ 

-=\ogN, ij, = -logg (144) 

we conclude that 

where Z{ij., k) is defined by eq. (59)-(60) and 

Z{g,N) = J d<p exp (-^tr^^ + ^tr,^^^ . (146) 

The integral (146) is of course divergent and should just be regarded as a 
shorthand for the formal power series in the coupling constant g. 

It is straightforward to generalize the preceding arguments to the case 
of general unrestricted triangulations, where arbitrary polygons are allowed. 
Instead of one coupling constant g we have a set gi, §2, gs, ■ ■ ■, but we will still 
use the notation g or gi. In this case cq. (146) is replaced by 



Z(5i,7V))= / ^e-^*'-^^'^), (147) 



where the potential V, which depends on all the coupling constants gi, is given 
by 

In eq. (148) we have scaled (p ^/N(j> for later convenience. It is not difficult 
to check that the obvious generalization of eq. (145) holds and the weight of 
a triangulation T is given by 

c^'m(^)J{g^^^^\ 
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where Nj{T) is the number of j'-gons in T. eq. (147) is of course a representa- 
tion of a formal power scries which is obtained by expanding the exponential 
of the non-quadratic terms as a power series in the coupling constants and 
then performing the Gaussian integrations term by term. 

Differentiating \ogZ{gi, N) with respect to the coupling constants g^, one 
obtains the expectation values of products of traces of powers of (f). These 
expectation values have a straightforward interpretation in terms of triangu- 
lations. Denoting the expectation with respect to the measure 

by (•) wc sec, for example, that (^N~^tr cf)"'^ is given by the sum over all 
connected triangulations of arbitrary genus whose boundary is an n-gon with 
one marked link. Similarly, 

^ (tr 0"tr 0™) - ^ (tr 0") (tr 0™) (149) 

is given by the sum over all connected triangulations whose boundary consists 
of two components with n and m links. More generally, the relation to the 
combinatorial problem is given by 

ki,...,kt ^1 '"^b 

where the subscript conn indicates the connected part of the expectation (•). 
One can rewrite eq. (150) as 

w{gi,zi,...,Zb) = N'^-^tv^— ■■■ tr^—) . (151) 

The one-loop function w{gi, z) is related to the density p{\) of eigenvalues 
of (j) defined by 

p{\) = l^5{\-\)Y (152) 

where Aj, i = 1,. . . ,N denote the A'' eigenvalues of the matrix (j). With this 
definition we have 

-i(tr0")= / rfAp(A)A", n>0. (153) 

■'■^ J—00 

Hence, 

w{g^,z)^ r dX-^. (154) 
J-00 z — A 

In the limit N 00 the support of p is confined to a finite interval [c_ , c+] 
on the real axis. In this case 'w{z) will be an analytic function in the complex 
plane, except for a cut along the interval [c_ , c+] . Note that p(A) is determined 
from w{z) by 

27rip(A) = lim (w;(A - is) - w{X + ie)) . (155) 

e— >0 
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6.1 The loop equations 

A standard method in quantum field theory is to derive identities by a change 
of variables in functional integrals. Here we apply this method to the matrix 
models and explore the invariance of the matrix integral (147) under infinites- 
imal field redefinitions of the form 

^^0 + e?i", (156) 

where e is an infinitesimal parameter. One can show that to first order in e 
the measure d<p defined by eq. (141) transforms as 

#^ # j 1 + e^tr/ trc^i"-*^ I . (157) 
V fe=o / 

The action transforms according to 

tvV{<j))^tYV{(l))+etr<i)''V'{(l)) (158) 

to first order in e. We can use these formulas to study the transformation of 
the measure under more general field redefinitions of the form 

fe=0 ^ 

This field redefinition only makes sense if z is on the real axis outside the 
support p. In the limit N ^ oo this is possible for z outside the interval 
[c_, c+]. Under the field redefinitions (159) the transformations of the measure 
and the action are given by 

l + £tr^- tr^— - I , (160) 
Z — (f) z — (j) J 

tr V{(l>)^trV {(!))+ e (—^V\(t>)\ . (161) 




The integral (147) is of course invariant under this change of the integration 
variables. By use of eqs. (160) and (161) we obtain the identity 

. A^tr (^-L_v'{c!>)^ I e-^*'-^^^) = 0. (162) 

The contribution to the integral coming from the first term in {•} in eq. (162) 
is, by definition, 

N'^w'^{z)+w{z,z). (163) 

The contribution from the second term inside {•} in eq. (162) can be written 
as an integral over the one- loop function as follows: 
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N\ z-(f>/ J '^^'z-X Jc2iTiz-uj ^" ^ ' 

where the second equality follows from eq. (155). The curve C encloses the 
support of p but not z. It is essential for the existence of C that p have compact 
support. We can then write cq. (162) in the form 

f du! V'(u)) , , 9, , 1 , , ,,„„s 

where z is outside the interval [c_, c+] on the real axis. Since both sides of eq. 
(165) can be analytically continued to C \ [c_, c+] the equation holds in this 
domain. 

Wc recognize eq. (165) as the loop equation already derived by combina- 
torial means, except that wc here have an additional term involving w{z, z) 
and with a coefficient 1/iV^. In fact. (165) is the starting point for a 1/7V^- 
expansion, i.e. a higher-genus expansion. To leading order in 1/iV^ we have 
(as already derived) 

w{z) = wo{z) = \ (y'{z) - M{zW{z - c+){z - c_)) (166) 
and from eq. (155) the corresponding eigenvalue density is 



p(A) = — M(A)V(c+-A)(A-c_). (167) 

Let us now discuss how the new scaling limit is described in the matrix 
formalism. Hcrmitian matrix models arc often analyzed in terms of the dy- 
namics of their eigenvalues. Since the action in (147) is invariant under the 
transformation (p U(t>W , with U € U{N) a unitary N x A^-matrix, one can 
integrate out the "angular" degrees of freedom. What is left is an integration 
over the eigenvalues Aj of ^ only, 

N 

Z{g) oc l[dXi e-^^.-^(^^) J] ' ^H', (168) 

•' i=l k<l 

where the last factor, the Vandermonde determinant, comes from integrating 
over the angular variables, and where 

tvV{ct>)=J2n>^i)- (169) 

i=l 

Naively one might expect that the large- A'^ limit is dominated by a saddle- 
point with V'{X) = 0. However, this is not the case since the Vandermonde 
determinant in (168) contributes in the large- limit. The cut which appears 
in ^^{z) is a direct result of the presence of the Vandermonde determinant. In 
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this way one can say that the dynamics of the eigenvalues is "non-classical" , 
deviating from V'{X) = 0, the size of the cut being a measure of this non- 
classicality. To get to the generalized CDT model we introduced a new cou- 
pling constant Qs in the matrix model by substituting 

y(^) ^ 1 v{cl>), (170) 
9s 

and considered the limit gs — > 0. As we have seen the coupling constant Qs 

controls and reduces the size of the cut and thus brings the system closer to 
a "classical" behaviour. Thus the quantum fluctuations are reduced in the 
generalized CDT models. 

Let us now consider the matrix potential (116), which formed the starting 
point of our new scaling analysis. We are still free to perform a change of 
variables. Inspired by relations (123)-(126), let us transform to new "CDT"- 
variables 

(j)^ zi + a$ + 0{a^), (171) 
at the same time re-expressing g as 

g = g{l- a2yl,dt + 0(a^)), (172) 

following eq. (123). Substituting the variable change into the matrix potential, 
and discarding a ^-independent constant term, one obtains 

A ,.(6 _ 1(53 

Vi<t>) = m ^ '^^^^^^ ' (173) 
in the limit a 0, from which it follows that 

Z{g,g,) = a^'z(.l,dt, G,), ZiA.^t, G,) = J ^-mrn^)_ (^74) 

The disk amplitude for the potential V{<P) is precisely W{Zcdt, ^cdt, Gs), and 
since by definition 

111 , , 

(175) 



z - (t> a Zcdt - ^ ' 

the first equal sign in eq. (129) follows straightforwardly from the simple 
algebraic equation (173). We conclude that the continuum generalized CDT- 
theory is described by the matrix model with potential V{$). 

The conclusion is that to leading order in N, the combinatorial method 
which works with a regularized lattice theory with a geometric interpreta- 
tion and an explicit cut-off a has a matrix model representation, and even 
in the continuum limit where a there exists a matrix model representa- 
tion of the theory. Once this is proved, one can actually "derive" a number 
of the results known for the generalized CDT model from the matrix model 
by the formal manipulations of eqs. (173)- (174), where a appears merely as 
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a parameter without any obvious geometric interpretation as a cut-off. Also, 
once the matrix model equivalence is established, it is clear that one has au- 
tomatically a higher-genus expansion available. This higher-genus expansion 
was first established working with the combinatorial, generalized CDT theory, 
and based on purely geometric arguments of splitting and joining of space as 
a function of proper time, as we (partly) described above. Such a theory could 
be turned into a kind of string field theory, imitating the work of Kawai and 
collaborators on Euclidean quantum gravity [34]. Like in the Euclidean case 
of Kawai et al., the formulation of the CDT string field theory used entirely 
a continuum notation (i.e. to cut-off a was already taken to zero). However, 
contrary to the situation for the Euclidean string field theory, we now have a 
matrix model representation even in the continuum, and many of the string 
field theory results derived in [17] follow easily from the loop equations (165) 
for the potential (173). 



6.2 Summation over all genera in the CDT matrix model 

One remarkable application of the CDT matrix model representation is that 
we can find the Hartle-Hawking wave function summed over all genera. 

In [17] the matrix model given by eqs. (173) and (174) was related to 
the CDT Dyson- Schwinger equations by (i) introducing into the latter an 
expansion parameter a, which kept track of the genus of the two-dimensional 
spacctime, and (ii) identifying this parameter with l/N"^. where N is the size 
of the matrix in the matrix integral. The 1/iV-expansion of our matrix model 
therefore plays a role similar to the l/A''-expansion originally introduced by 
't Hooft: it reorganizes an asymptotic expansion in a coupling constant t 
{t = Gs/{^/ A^AtY^^ in our case) into convergent sub-summations in which 
the fcth summand appears with a coefficient N~'^''. In QCD applications, 
the physically relevant value is iV = 3, to which the leading-order terms in 
the large iV-expansion can under favourable circumstances give a reasonable 
approximation. 

As we will see, for the purposes of solving our string field-theoretic model 
non-perturbatively, an additional expansion in inverse powers of N (and thus 
an identification of the contributions at each particular genus) is neither es- 
sential nor does it provide any new insights. This means that we will consider 
the entire sum over topologies "in one go" , which simply amounts to setting 
N = 1, upon which the matrix integral (174) reduces to the ordinary integral^ 



Z{Gs,Acdt) = J dm exp ^^cdtm - ^ m^^ 



(178) 



Starting from a matrix integral for N x AT-matrices like (174), performing a formal 
expansion in (matrix) powers commutes with setting A^ = 1, as follows from the 
following property of expectation values of products of traces, which holds for 
any n = 1, 2, 3, . . . and any set of non-negative integers {rik}, k = 1, . . . , 2n, such 
that X]fe=i = 2n. For any particular choice of such numbers, consider 
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while the disk amplitude can be written as 

WMX) = . . / dm ^ ^. (179) 

Z{GB,Acdt)J X-m 

These integrals should be understood as formal power series in the dimension- 
less variable t = Gg/iVAcdt)'^^^ appearing in eq. (128). Any choice of an inte- 
gration contour which makes the integral well defined and reproduces the for- 
mal power series is a potential non-perturbative definition. However, difi^erent 
contours might produce different non-perturbative contributions (i.e. which 
cannot be expanded in powers of t), and there may even be non-perturbative 
contributions which are not captured by any choice of integration contour. 
As usual in such situations, additional physics input is needed to fix these 
contributions. 

To illustrate the point, let us start by evaluating the partition function 
given in (178). We have to decide on an integration path in the complex plane 
in order to define the integral. One possibility is to take a path along the 
negative axis and then along either the positive or the negative imaginary 
axis. The corresponding integrals are 

Z{g, A) = t'/'F±{t-^/'), F±(t-2/3) = 277 e±^^/'^Ai(t-2/3e±2-/3)^ 

(180) 

where Ai denotes the Airy function. Both F± have the same asymptotic ex- 
pansion in t, with positive coefficients. Had we chosen the integration path en- 
tirely along the imaginary axis we would have obtained (27ri times) Ai(t~^/^), 
but this has an asymptotic expansion in t with coefiicients of oscillating sign, 



/ dM e-- ^'^ n^:! (tr /n) 



J2 ^mN"", (176) 



ra — — n 



where the last equation defines the numbers u)m as coefiicients in the power ex- 
pansion in N of the expectation value. Now, we have that 



Wm = (2n-1)!! (177) 



m=—'n 



independent of the choice of partition {uk}- The number (2n— 1)!! simply counts 
the "Wick contractions" of a;^" which we could have obtained directly as the 
expectation value {x^"), evaluated with a one-dimensional Gaussian measure. 
In the model at hand, wc will calculate sums of the form X]m=-n directly, 
since we are summing over all genera without introducing an additional coupling 
constant for the genus expansion. In other words, the dimensionless coupling 
constant t in this case already contains the information about the splitting and 
joining of the surfaces, and the coefficient of contains contributions from 2d 
geometries whose genus ranges between and [k/2]. Wc cannot disentangle these 
contributions further unless we introduce A'^ as an extra parameter. 
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which is at odds with its expHcit power expansion in t. We have (using the 
standard notation of Airy functions) 

F±(^) = 7r(Bi(2) ±iAi(^)), (181) 

from which one deduces immediately that the functions F±(t~'^^^) are not 
real. However, since Bi{t^^/-^) grows like est for small t while Ai(t~^/^) falls 
off like e~3t , their imaginary parts are exponentially small in l/t compared 
to the real part, and therefore do not contribute to the asymptotic expansion 
in t. An obvious way to define a partition function which is real and shares 
the same asymptotic expansion is by symmetrization, 

i(F+ + f_)^^Bi. (182) 

The situation parallels the one encountered in the double-scaling limit of the 
"old" matrix model but is less complicated. 

Presently, let us collectively denote by F{z) any of the functions F±{z) or 
7rBi(2;), leading to the tentative identification 

Z{G,, vlcdt) = t'/^ F{t-^'^) , F"{z) = zF{z), (183) 

where we have included the differential equation satisfied by the Airy functions 
for later reference. Assuming X > 0, we can write 

1 f°° 

— = / dL exp [- {X -m)L]. (184) 

X -m Jq 

We can use this identity in eq. (179) to obtain the integral representation 

F(t-y^-t'/^^,L) 



W^cdt(^) = / dL e-^^ ^ y ^ ^. (185) 

Jo F f-^/A 



From the explicit expression of the Laplace transform we can now read off the 
Hartle-Hawking amplitude as function of the boundary length L, 

Before turning to a discussion of the non-perturbative expression for Wcdt(-^) 
we have just derived, let us remark that the asymptotic expansion in t of course 
agrees with that obtained by recursively solving the CDT Dyson- Schwinger 
equations. Using the standard asymptotic expansion of the Airy function one 
obtains 



W^,,,(L)=e-v^^e*''(*'^^^) ^k^oCkf'O^^tVA^.L) ^^^^^ 

Z^k=0 ^k ^ 
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where the coefRcients Ck arc given by cq — 1, Ck ^ ^ (j)'' (^) j, fc > 0. In 

(187), we have rearranged the exponential factors to exhibit the exponential 
fall-off in the length variable L, multiplied by a term containing the function 



2 



(1 - t^/A^tLf/^ - 1 + It^/A^tL 



(188) 



which has an expansion in positive powers of t. 

Wcdt(i) has the interpretation of the wave function of the spatial uni- 
verse according to the hypothesis of Hartle and Hawking. L € [0, oo] and the 
probability of finding a spatial universe with length between L and L -|- dL is 

P(L) = (189) 

since the integration measure is dL/L. Thus the probability is not normal- 
izable in a conventional way and peaked at L = since Wcdt{L = 0) = 1. 
However, for each term in the asymptotic expansion (187) we obtain a finite 
value (L) ~ 1/V^cdt as one would naturally expect. Since termination of the 
series in (187) at a finite k also implies restricting the 2d spacetime to have 
a finite genus we can say that as long as we restrict spacetime to have finite 
genus we have (L) ~ 1/V^cdt- However, if we allow spacetimes of arbitrarily 
large genus to appear, i.e. if topology fluctuations are unconstrained (that 
means at most suppressed by a coupling constant, but no upper limit on the 
genus imposed by hand), a remarkable change appears: {L) = oo because the 
full non-perturbative Wcdt (L) does not fall off like e~^'^<='^' ^ but only as L~^/'^ 
(Note that Wcdt {L) is still integrable at infinity since the integration meastire 
is dL/L). This dramatic change in large L behaviour {W{L) also becomes 
oscillatory for large L, despite the fact that each term in the asymptotic ex- 
pansion (187) is positive) is clearly to be attributed to surfaces of arbitrarily 
large genus, i.e. it is a genuinely non-perturbative result. 



7 Discussion and perspectives 

The four-dimensional CDT model of quantum gravity is extremely simple. 
It is the path integral over the class of causal geometries with a global time 
foliation. In order to perform the summation explicitly, we introduce a grid 
of piecewise linear geometries, much in the same way as when defining the 
path integral in quantum mechanics. Next, we rotate each of these geometries 
to Euclidean signature and use as bare action the Einstein-Hilbert action* in 
Regge form. That is all. 



Of course, the full, effective action, including measure contributions, will contain 
all higher-derivative terms. 
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The resulting superposition exhibits a nontrivial scaling behaviour as func- 
tion of the four-volume, and we observe the appearance of a well-defined av- 
erage geometry, that of de Sitter space, the maximally symmetric solution 
to the classical Einstein equations in the presence of a positive cosmological 
constant. We are definitely in a quantum regime, since the fluctuations of 
the three-volume around de Sitter space are sizable, as can be seen in Fig. 7. 
Both the average geometry and the quantum fluctuations are well described 
in terms of the minisuperspace action (25). A key feature to appreciate is 
that, unlike in standard (quantum-) cosmological treatments, this description 
is the outcome of a non-perturbative evaluation of the full path integral, with 
everything but the scale factor (equivalently, V^{t)) summed over. Measuring 
the correlations of the quantum fluctuations in the computer simulations for a 
particular choice of bare coupling constants enabled us to determine the con- 
tinuum gravitational coupling constant G as G ~ 0.42a^, thereby introducing 
an absolute physical length scale into the dimensionless lattice setting. Within 
measuring accuracy, our de Sitter universes (with volumes lying in the range 
of 6.000-47.000 £pi) are seen to behave perfectly semi-classically with regard 
to their large-scale properties. 

We have also indicated how we may be able to penetrate into the sub- 
Planckian regime by suitably changing the bare coupling constants. By "sub- 
Planckian regime" we mean that the lattice spacing a is (much) smaller than 
the Planck length. While we have not yet analyzed this region in detail, we 
expect to eventually observe a breakdown of the semi-classical approximation. 
This will hopefully allow us to make contact with continuum attempts to 
define a theory of quantum gravity based on quantum field theory. One such 
attempt has been described in the introduction and is based on the concept of 
asymptotic safety. It uses renormalization group techniques in the continuum 
to study scaling violations in quantum gravity around a UV fixed point [2]. 
Other recent continuum field theoretical models of quantum gravity which are 
not in disagreement with our data are the so-called Lifshitz gravity model [3] 
and the so-called scale-invariant gravity model [5, 4]. In principle it is only 
a question of computer power to decide if any of the models agree with our 
CDT model of quantum gravity. 

On the basis of these results two major issues suggest themselves for further 
research. First, we need to establish the relation of our effective gravitational 
coupling constant G with a more conventional gravitational coupling constant, 
defined directly in terms of coupling matter to gravity. In the present work, 
we have defined G as the coiipling constant in front of the effective action, but 
it would be desirable to verify directly that a gravitational coupling defined 
via the coupling to matter agrees with our G. In principle it is easy to couple 
matter to our model, but it is less straightforward to define in a simple way a 
set-up for extracting the semi-classical effect of gravity on the matter sector. 
Attempts in this direction were already undertaken in the "old" Euclidean 
approach [29, 30], and it is possible that similar ideas can be used in CDT 
quantum gravity. 
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The second issue concerns the precise nature of the "continuum limit". 
Recall our discussion in the Introduction about this in a conventional lattice- 
theoretic setting. The continuum limit is usually linked to a divergent correla- 
tion length at a critical point. It is unclear whether such a scenario is realized 
in our case. In general, it is rather unclear how one could define at all the con- 
cept of a divergent length related to correlators in quantum gravity, since one 
is integrating over all geometries, and it is the geometries which dynamically 
give rise to the notion of "length" . 

This has been studied in detail in two-dimensional (Euclidean) quantum 
gravity coupled to matter with central charge c < 1 [31]. It led to the conclu- 
sion that one could associate the critical behaviour of the matter fields (i.e. 
approaching the critical point of the Ising model) with a divergent correla- 
tion length, although the matter correlators themselves had to be defined as 
non-local objects due to the requirement of diffeomorphism invariancc. On 
the other hand, the two-dimensional studies do not give us a clue of how 
to treat the gravitational sector itself, since they do not possess gravitational 
field-theoretic degrees of freedom. As we have seen the two-dimensional lattice 
models can be solved analytically and the only fine-tuning needed to approach 
the continuum limit is an additive renormalization of the cosmological con- 
stant. Thus, fixing the two-dimensional spacetime volume N2 (the number 
of triangles), such that the cosmological constant plays no role, there are no 
further coupling constants to adjust and the continuum limit is automatically 
obtained by the assignment V2 = A^2a^ and taking A^2 ^ 00. This situation 
can also occur in special circumstances in ordinary lattice field theory. A term 



(or a higher-dimensional generalization) will also go to the continuum free 
field theory simply by increasing the lattice size and using the identification 
Vd = L'^a'^ {L denoting the linear size of the lattice in lattice units), the 
higher-derivative term being sub-dominant in the limit. It is not obvious that 
in quantum gravity one can obtain a continuum quantum field theory with- 
out fine-tuning in a similar way, because the action in this case is multiplied 
by a dimensionful coupling constant. Nevertheless, it is certainly remarkable 
that the infrared limit of our effective action apparently reproduces - within 
the cosmological setting - the Einstein-Hilbert action, which is the unique 
diffeomorphism-invariant generalization of the ordinary kinetic term, contain- 
ing at most second derivatives of the metric. A major question is whether 
and how far our theory can be pushed towards an ultraviolet limit. We have 
indicated how to obtain such a limit by varying the bare coupling constants 
of the theory, but the investigation of the limit a ^ with fixed G has only 
just begun and other scenarios than a conventional UV fixed point might be 
possible. One scenario, which has often been discussed as a possibility, but 
which is still missing an explicit implementation is the following: when one 
approaches sub-Planckian scales the theory effectively becomes a topological 



like 




(190) 
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quantum field theory where the metric plays no role. Also in our very explicit 
implementation of a quantum gravity model it is unclear how such a scenario 
would look. 
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